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DECLARATION OF AVI ASHKENAZL Ph,D UNDER 37 C.F.R. $ 1.132 

I, Avi Ashkenazi, Ph.D. declare and say as follows: - 

1 . I am Director and Staff Scientist at the Molecular Oncology Department of 
Genentech, Inc., South San Francisco, CA 94080. 

2. Ijoined Genentech in 1988 as a postdoctoral fellow. Since then, I have 
investigated a variety of cellular signal transduction mechanisms, including apoptosis, and have 
developed technologies to modulate such mechanisms as a means of therapeutic intervention in 
cancer and autoimmune disease. I am currently involved in the investigation of a series of 
secreted proteins over-expressed in tumors, with the aim to identify useful targets for the 
development of therapeutic antibodies for cancer treatment. 

3. My scientific Curriculum Vitae, including my Ust of publications, is attached to 
and forms part of this Declaration (Exhibit A). 

4. Gene amplification is a process in which chromosomes undergo changes to 
contain multiple copies of certain genes that normally exist as a single copy, and is an important 
factor in the pathophysiology of cancer. Amplification of certain genes (e.g., Myc or Her2/Neu) 



gives cancer cells a growth or survival advantage relative to normal cells, and might also provide 
a mechanism of tumor cell resistance to chemotherapy or radiotherapy. 

5. Ifgene amplification results in over-expression ofthemRNA and the 

corresponding gene product/then it identifies that gene product as a promising target for cancer 
therapy, for example by the therapeutic antibody approach. Even in the absence of over- 
expression of the gene product, amplification of a cancer marker gene - as detected, for example, 
by the reverse transcriptase TaqMan® PGR or the fluorescence in situ hybridization (FISH) 
assays -is usefiil in the diagnosis or classification of cancer, or in predicting or monitoring the 
efficacy of cancer therapy. An increase in gene copy number can result not only from 
intrachromosomal changes but also from chromosomal aneuploidy. It is important to understand 
that detection of gene ampUfication can be used for cancer diagnosis even if the determination 
includes ineasurement of chromosomal aneuploidy. Indeed, as long as a significant difference 
relative to normal tissue is detected, it is irrelevant if the signal originates from an increase in the 
number of gene copies per chromosome and/or an abnormal number of chromosomes. 

6. I understand that according to the Patent Office, absent data demonstrating that 
the increased copy number of a gene m certain types of cancer leads to increased expression of 
its product, gene amplification data are insufficient to provide substantial utility or well 
established utiUty for the gene product (the encoded polypeptide), or an antibody specifically 
binding the encoded polypeptide. However, even when ampUfication of a cancer marker gene 
does not result in significant over-expression of the corresponding gene product, this very 
absence of gene product over-expression still provides significant mformation for cancer 
diagnosis and freatment. Thus, if over-expression of the gene product does not parallel gene 
ampUfication in certain tumor types but does so in others, then parallel monitoring of gene 
ampUfication and gene product over-expression enables more accurate tumor classification and 
hence better determination of suitable therapy. In addition, absence of over-expression is crucial 
information for the practicing cUnician. If a gene is amplified but the corresponding gene 

. product is not over-expressed, the cUnician accordingly wiU decide riot to treat a patient with 
agents that target that gene product. 

7 . I hereby declare that all statements made herein of my own knowledge are true 
and that aU statements made on information or beUef are beUeved to be true, and fiirther that 
these statements were, made with the knowledge that willfiil false statements and the like so 



made are punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the 
United States Code and that such willful statements may jeopardize the vahdity of the 
appUcation or any patent issued diereon. 
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1456 Tarrytown Street, San Mateo, CA 94402 
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Education: 

1983: 
1986: 



B.S. in Biochemiistry, with honors, Hebrew University, Israel 
Ph^D. in Biochemistry, Hebrew University, Israel 
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1985- 1986: 

1986- 1988: 

1988- 1989: 

1989- 1993: 
1994-1996: 

1996- 1997: 

1997- 1990: 
1999-2002: 
2002-present: 



Teaching assistant, undergraduate level course in Biochemistry 
Teaching assistant, graduate level course on Signal Transduction 
Postdoctoral fellow, Hormone Research Dept., UCSF, and 
Developmental Biology Dept., Genentech, hic, with J. Ramachandran 
Postdoctoral fellow. Molecular Biology Dept., Genentech, Inc., 
with D. Capon 

Scientist, Molecular Biology Dept., Genentech, lac. 
Senior Scientist, Molecular Oncology Dept., Genentech, Inc. 
Senior Scientist and Interim director, Molecular Oncology Dept., 
Genentech, Inc. 

Senior Scientist and preclinical project team leader, Genentech, Inc. 

Staff Scientist in Molecular Oncology, Genentech, Inc. 

Staff Scientist and Director in Molecular Oncology, Genentech, Inc. 
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1988: 



First prize. The Boehringer Ingelheim Award 



1 



Editorial: 

Editorial Board Member: Current Biology 
Associate Editor, Clinical Cancer Research. 
Associate Editor, Cancer Biology and Therapy. 

Refereed papers: 

1 . Gertler, A., Ashkenazi. A., and Madar, Z. Binding sites for human growth 
hormone and ovine and bovine prolactins in the mammary gland and liver of the 
lactating cow. Mol. Cell. Endocrinol. 34, 5U57 (1984). 

2. Gertler, A., Shamay, A., Cohen, N., Ashkenazi, A., Friesen, H., Levanon, A., 
Gorecki, M., Aviv, H., Hadari, D., and Vogel, T. hihibition of lactogenic 
activities of ovine prolactm and human growth hormone (hGH) by a novel form of 

a modified recombinant hGH. Endocrinology 118, 720-726 (1986). 

3. Ashkenazi. A., Madar, Z., and Gertler, A. Partial purification and characterization 
of bovine mammary gland prolactin receptor. Mol. Cell. Endocrinol. 50, 79-87 

- (1987). 

4. Ashkenazi. A., Pines, M., and Gertler, A. Down-regulation of lactogenic 
hormone receptors in Nb2 lymphoma cells by cholera toxin. Biochemistry 
Intematl 14, 1065-1072 (1987). 

5. Ashkenazi. A.. Cohen, R., and Gertler, A. Characterization of lactogen receptors 
in lactogenic hormone-dependent and independent Nb2 lymphoma cell lines. 
F£:5SIe«. 210, 51-55 (1987). 

6. Ashkenazi. A., Vogel, T., Barash, I., Hadari, D., Levanon, A.. Gorecki, M., and 
Gertler, A. Comparative study on in vitro and in vivo modulation of lactogemc 
and somatotropic receptors by native human growth hormone and its modified 
recombinant analog. i:«i/ocn>io/ogy 121, 414-419 (1987). 

7. Peralta, E., Winslow, J., Peterson, G., Smith, D., Ashkenazi, A., Ramachandran, 
J., Schimerlik, M., and Capon, D. Primary structure and biochemical properties 
of an M2 muscarinic receptor. Science 236, 600-605 (1987). 

8. Peralta, E. Ashkenazi. A.. Winslow, J., Smith, D., Ramachandran, J., and Capon, 
D. J. Distincnt primary structures, Ugand-binding properties and tissue-specific 

expression of four human muscarinic acetylcholine receptors. EMBO J. 6, 3923- 

3929 (1987). 

9. Ashkenazi. A. . Winslow, J., Peralta, E., Peterson, G., Schimerlik, M., Capon, D., 
and Ramachandran, J. An M2 muscarinic receptor subtype coupled to both 
adenylyl cyclase and phosphoinositide turnover. Science 238, 672-675 (1987). 



10. Pines, M., Ashkenazi, A., Cohen-Chapnik, N., Binder, L., and Gertler, A. 
Inhibition of the proliferation of Nb2 lymphoma cells by femtomolar 
concentrations of cholera toxin and partial reversal of the effect by 1 2-o- 
tetradecanoyl-phorbol-13-acetate. J. Cell. Biochem. 37, 1 19-129 (1988). 

1 1 . Peralta, E. Ashkenazi. A.. Winslow, J. Ramachandran, J., and Capon, D. 
Differential regulation of PI hydrolysis and adehylyl cyclase by muscarinic 
receptor subtypes. T/a^wre 334, 434-437 (1988). 

12. Ashkenazi.. A. Peralta, E., Winslow, J., Ramachandran, J ., and Capon, D. 
Functionally distinct G proteins couple different receptors to PI hydrolysis in the 
same cell. Cell 56, 487-493 (1989). 

13. Ashkenazi. A., Ramachandran, J., and Capon, D. Acetylcholine analogue 
stimulates DNA synthesis in brain-derived ceUs via specific muscarinic 
acetylcholine receptor subtypes. Nature 340, 146-150 (1989). 

14. Lammare, D., Ashkenazi. A., Fleury, S., Smith, D., Sekaly, R., and Capon, D. 
The MHC-binding and gpl20-binding domains of CD4 are distinct and separable. 
Science 245, 743-745 (1989). 

15. Ashkenazi.. A.. Presta, L., Marsters, S., Camerato, T., Rosenthal, K., Fendly, B., 
and Capon, D. Mapping the CD4 bmding site for human immunodefficiency 
virus type 1 by alanine-scanning mutiagenesis. Prop. Natl. Acad. Sci. USA. 87, 
7150-7154(1990). 

16. Chamow, S., Peers, D., Byra, R., Mulkeirin, M., Harris, R., Wang, W.,Bjorkman, 
P., Capon, D., and Ashkenazi. A. Enzymatic cleavage of a CD4 immunoadhesin 
generates crystallizable, biologically active Fd-like fragments. Biochemistry 29, 
9885-9891 (199.0). 

17. Ashkenazi. A.. Smith, D., Marsters, S., Riddle, L., Gregory, T., Ho, D., and 
Capon, D. Resistance of primary isolates of human immunodefficiency virus type 
1 to soluble CD4 is independent of CD4-rgp 120 binding affinity. Proc. Natl. 
Acad. Sci. USA. 88, 7056-7060 (1991). 

18. Ashkenazi. A.. Marsters, S., Capon, D., Chamow, S., Figari., I., Pennica, D., 
Goeddel., D., Palladino, M., and Smith, D. Protection against endotoxic shock by 
a tumor necrosis factor receptor immunoadhesin. Proc. Natl. Acad. Sci. USA. 88, 
10535-10539(1991). 

19. Moore, J., McKeating, J., Huang, Y., Miken^LA,., and Ho, D. Virions of 
primary HIV-1 isolates resistant to sCD4 neutralization differ in sCD4 affinity and 
glycoprotein gpl20 retention from sCD4-sensitive isolates. J. Virol. 66, 235-243 
(1992). 



3 



20. Jin, H., Oksenberg, D., Ashkenazi. A.. Peroutka, S., Duncan, A., Rozraahel., R., 
Yang, Y., Mengod, G., Palacios, J., and O'Dowd, B. Characterization of the 
human 5-hydroxytryptamineiB receptor. J. Biol. Chem. 261, 5735-5738 (1992). 

21. Marsters, A., Frutkin, A., Sunpson, N., Fendly, B. and Ashkenazi. A. 
Identification of cysteine-rich domains of the type 1 tumor necrosis receptor 
involved in hgand binding. Biol. Chem. 267, 5747-5750 (1992). 

22. Chamow, S., Kogan, T., Peers, D., Hastings, R., Bym, R., and Ashkenazi. A. 
Conjugation of sCD4 without loss of biological activity via a novel carbohydrate- 
directed cross-linking reagent. J. Biol. Chem. 267, 15916-15922 (1992). 

23. Oksenberg, D., Marsters, A., O'Dowd, B., Jin, H., Havlik, S., Peroutka, S., and 
Ashkenazi. A. A single amino-acid difference confers major pharmacologic 
variation between human and rodent 5-HTiB receptors. Nature 360, 161-163 

(1992) . 

24. Haak-Frendscho, M., Marsters, S., Chamow, S., Peers, D., Sunpson, N., and 
Ashkenazi. A. Inhibition of interferon y by an interferon y receptor 
immunoadhesin. Immunology 79, 594-599 (1993). 

25. Penica, D., Lam, V., Weber, R., Kohr, W., Basa, L., Spelhnan, M., Ashkenazi. 
Shire, S ., and Goeddel, D. Biochemical characterization of the extracellular 
domain of the 75-kd tumor necrosis factor receptor. Biochemistry 32, 3131-3138. 

(1993) . 

26. Barfod, L., Zheng, Y., Kuang, W., Hart, M., Evans; T., Cerione, R., and 
Ashkenazi. A. Cloning and expression of a human CDC42 GTPase Activating 
Protein reveals a functional SH3-binding domain. J. Biol. Chem. 268, 26059- 
26062(1993): 

27. Chamow, S., Zhang, D., Tan, X., Mhtire, S., Marsters, S., Peers, D., Bym, R., 
Ashkenazi. A., and Yunghans, R. A humanized bispecific immunoadhesin- 
antibody that retargets CD3+ effectors to kill HIV- 1 -infected cells. J. Immunol 
153,4268-4280(1994). 

28. Means, R., Rrahtz, S., Luna, J., Marsters, S., and Ashkenazi. A. Inhibition of 
murine erythroid colony formation in vitro by iterferon y and correction by 
interferon y receptor immunoadhesin. Blood 83, 911-915 (1994). 

29. Haak-Frendscho, M., Marsters, S., Mordenti, J., Gillet, N., Chen, S., 
an dAshkenazi. A. hihibition of TNF by a TNF receptor inmiunoadhesin: 
comparison with an anti-TNF mAb. J. Immunol. 152, 1347-1353 (1994). 
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30. Chamow, S., Kogan, T., Venuti, M., Gadek, T., Peers, D., Mordenti, J., Shak, S., 
and Ashkenazi. A. Modification of CD4 immunoadhesin with monomethoxy- 
PEG aldehyde via reductive alkilation. Bioconj. Chem. 5, 133-140 (1994). 

31. Jin, H., Yang, R., Marsters, S., Bunting, S., Wurm, F., Chamow, S., and 
Ashkenazi, A. Protection against rat endptoxic shock by p55 tumor necrosis factor 
(TNF) receptor immunoadhesin: comparison to anti-TNF monoclonal antibody. J. 
Infect. Diseases 170, 1323-1326 (1994). 

32. Beck, J., Marsters, S., Harris, R., Ashkenazi. A., and Chamow, S. Generation of 
soluble interleukin-1 receptor from an immunoadhesin by specific cleavage. Mol. 
/7n/«M«o/. 31,1335-1344 (1994). 

33. Pitti, B., Marsters, M., Haak-Frendscho, M., Osaka, G., Mordenti, J., Chamow, S., 
and Ashkenazi, A. Molecular and biological properties of an mterleukin- 1 
receptor immunoadhesin. Mol. Immunol. 31, 1345-1351 (1994). 

34. Oksenberg, D., Havlik, S., Peroutka, S., and Ashkenazi. A. The third intracellular 
loop of the 5-HT2 receptor specifies effector coupling. /. Neurochem. 64, 1440- 
1447(1995). 

35. Bach, E., Szabo, S., Dighe, A., Ashkenazi. A.. Aguet, M., Murphy, K., and 
Schreiber, R. Ligand-induced autoregulation of IFN-y receptor P chain expression 
in T helper cell subsets. 5cie«ce 270, 1215-1218 (1995). 

36. Jin, H., Yang, R., Marsters, S., Ashkenazi. A. . Bunting, S., Marra, M., Scott, R., 
and Baker, J. Protection against endotoxic shock by bactericidal/permeability- 
increasing protein in rats. J. Clin. Invest. 95, 1947-1952 (1995). 

37. Marsters, S., Penica, D., Bach, E., Schreiber, R., and Ashkenazi. A. Interferon y 
signals via a high-affinity multisubunit receptor complex that contains two types 
of polypeptide chain: Proc. Natl. Acad. Sci. USA. 92, 5401-5405 (1995). 

38. Van Zee, K., Moldawer, L., Oldenburg, H., Thompson, W., Stackpole, S., 
Montegut, W., Rogy, M., Meschter, C, Gallati, H., Schiller, C, Richter, W., 
Loetcher, H., Ashkenazi. A .. Chamow, S., Wurm, F., Calvano, S., Lowry, S., and 
Lesslauer, W. Protection against lethal E. coU bacteremia in baboons by 
pretreatment with a 55-kDa TNF receptor-Ig fusion protein, Ro45-2081. J: 
Immunol. 156,2221-2230 (1996). 

39. Pitti, R., Marsteris, S., Ruppert, S., Donahue, C, Moore, A., and Ashkenazi. A. 
Induction of apoptosis by Apo-2 Ligand, a new member of the tumor necrosis 
factor cytokine family. /. Biol. Chem. 271, 12687-12690 (1996). 
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40. Marsters, S., Pitti, R., Donahue, C, Rupert, S., Bauer, K., and Ashkenazi. A. 
Activation of apoptosis by Apo-2 ligand is independent of FADD but blocked by 
CimA. Cwrr. 6, 1669-1676 (1996). 

41. . Marsters, S., Skubatch, M., Gray, C, and Ashkenazi. A . Herpesvirus entry 

mediator, a novel member of the tumor necrosis factor receptor family, activates 
the NF-kB and AP-1 transcription factors. Biol. Chem. 272, 14029-14032 
(1997). 

42. Sheridan, J., Marsters, S., Pitti, R., Gumey, A., Skubatch, M., Baldwin, D., 
Ramakrishnan, L., Gray, C, Baker, K., Wood, W.I., Goddard, A., Godowski, P., and 
Ashkenazi. A. Control of TRAIL-induced apoptosis by a family of signaling and 
decoy receptors, ^ciewce 277, 818-821 (1997). 

43. Marsters, S., Sheridan, J., Pitti, R., Gumey, A., Skubatch, M., Balswin, D., Huang, A., 
Yuan, J., Goddard, A., Godowski, P., and Ashkenazi. A. A novel receptor for 
Apo2L/TRAD:. contains a truncated death domain. Curr. Biol. 7, 1 003- 1 006 ( 1 997). 

44. Marsters, A., Sheridan, J., Pitti, R., Brush, J., Goddard, A., and Ashkenazi. A. 
Identification of a Ugand for the death-domain-contaming receptor Apo3. Curr. Biol. 
8,525-528 (1998). 

45. Rieger, J., Naumann, U., Glaser, T., Ashkenazi. A ., and Weller, M. Apo2 Ugand: 
a novel weapon against malignant glioma? FEBS Lett. 427, 124- 1 28 (1998). 

46. Pender, S., Fell, J., Chamow, S., Ashkenazi. A ., and MacDonald, T. A p55 TNF 
receptor immimoadhesin prevents T cell mediated intestinal injury by inhibiting 
matrix metalloproteinase production. /. Immunol. 160, 4098-4103 (1998). 

47. Pitti, R., Marsters, S., Lawrence, D., Roy, Kischkel, F., M., Dowd, P., Huang, A., 
Donahue, C, Sherwood, S., Baldwin, D., Godowski, P., Wood, W., Gumey, A., 
ffillan, K., Cohen, R., Goddard, A., Botstein, D., and MtaazL^ Genomic 
amplification of a decoy receptor for Fas ligand in lung and colon cancer. Nature 
396,699-703(1998). 

48. Mori, S., Marakami-Mori, K., Nakamura, S., Ashkenazi. A ., and Bonavida, B. 
Sensitization of AIDS Kaposi's sarcoma cells to Apo-2 Ugand-induced apoptosis 
by actinomycin D. y. /»JTO««o/. 162, 5616-5623 (1999). 

49. Gumey, A. Marsters, S., Huang, A., Pitti, R., Mark, M., Baldwin, D., Gray, A., 
Dowd, P., Bmsh, J., Heldens, S., Schow, P., Goddard, A., Wood, W., Baker, K., 
Godowski, P., and Ashkenazi. A.- Identification of a new member of the tumor 
necrosis factor family and its receptor, a human ortholog of mouse GFTR. Curr. 
Bio/. 9,215-218 (1999). 
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50. Ashkenazi. A ., Pai, R., Fong, s., Leung, S., Lawrence, D., Marsters, S., Blackie, 
C, Chang, L., McMurtrey, A., Hebert, A., DeForge, L., Khoumenis, I., Lewis, D., 
Harris, L., Bussiere, J., Koeppen, H., Shahrokh, Z., and Schwall, R. Safety and 
anti-tumor activity of recombinant soluble Apo2 ligand. J. Clin. Invest. 104, 155- 
162 (1999). 

51. Chuntharapai, A., Gibbs, V., Lu, J., Ow, A., Marsters, S., Ashkenazi, A., De Vos, 
A., Kim, K.J. Determination of residues involved in hgand binding and signal 
transmissiion in the human EFN-a receptor 2. J. Immunol. 163, 766-773 (1999). 

52. Johnsen, A.-C., Haux, J., Steinkjer, B., Nonstad, U., Egeberg, K., Sundan, A., 
Ashkenazi. A., and Espevik, T. Regulation of Apo2L/TRAIL expression in NK 
cells - involvement in NK cell-mediated cytotoxicity. Cytokine 1 1 , 664-672 
(1999). 

53. Roth, W., Iserimann, S., Naumann, U., Kugler, S., Bahr, M., Dichgans, 
Ashkenazi, A., and Weller, M. Eradication of mtracranial human maUgnant 
glioma xenografts by Apo2L/TRAIL. Biochem. Biophys. Res. Commun. 265, 479- 
483(1999). 

54. Hymowitz, S.G., Christinger, H.W., Fuh, G., Ultsch, M., O'Connell, M., Kelley, 
R.F., Ashkenazi. A. and de Vos, A.M. Triggering Cell Death: The Crystal 
Structure of Apo2L/TRAIL in a Complex with Death Receptor 5. Molec. Cell 4, 
563-571 (1999). 

55. Hymowitz, S.G., O'Connel, M.P., Utsch, M.H., Hurst, A., Totpal, K., Ashkenazi, 
Ay de Vos, A.M., Kelley, R.F. A unique zinc-binding site revealed by a high- 
resolution X-ray structure of homotrimeric Apo2L/TRAIL. Biochemistry 39, 633- 
640(2000). 

56. Zhou, Q., Fukushima, P., DeOraff, W., Mitchell, J.B., Stetler-Stevenson, M., 
Ashkenazi.'^A.. and Steeg, P.S. Radiation and the Apo2L/TRAIL apoptotic 
pathway preferentially inhibit the colonization of premaUgnant human breast 
cancer cells overexpressing cyclin D 1 . Cancer Res. 60, 26 1 1-26 1 5 (2000). 

57. Kischkel, F.C., Lawrence, D. A., Chuntharapai, A., Schow, P., Kim, J ., and 
Ashkenazi. A. Apo2L/TRAIL-dependent recruitment of endogenous FADD and 
Caspase-8 to death receptors 4 and 5. Immunity 12, 611-620 (2000). 

58. Yan, M., Marsters, S.A., Grewal, I.S., Wang, H., * Ashkenazi. A., and *Dixit, 
V.M. Identification of a receptor for BlyS demonstrates a crucial role in humoral 
immunity. Mz/wre /m/MMWo/. 1,37-41(2000). 
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59. Marsters, S.A., Yan, M., Pitti, R.M., Haas, P.E., Dixit, V.M., and Ashkenazi, A. 
Interaction of the TNF homologues BLyS and APRIL with the TNF receptor 
homologues BCMA and TACI. Cuir. Biol. 10, 785-788 (2000). 

60. Kischkei, F.C., and Ashkenazi. A . Combining enhanced metaboUc labeling with 
immunoblotting to detect interactions of endogenous cellular proteins. 
Biotechniques 29, 506-512 (2000). 

61. Lawrence, D., Shahrokh, Z., Marsters, S., Achilles, K., Shih, D. Mounhd, B., 
Hillan, K., Totpal, K. DeForge, L., Schow, P., Hooley, J., Sherwood, S., Pai, R., 
Leung, S., Khan, L., Gliniak, B., Bussiere, J., Smitti, C, Strom, S., Kelley, S., 
Fox, J., Thomas, D., and Ashkenazi. A. Differential hepatocyte toxicity of 
recombinant Apo2L/TRAIL versions. Nature Med. 7, 383-385 (2001). 

62. Chuntharapai, A., Dodge, K., Grimmer, K., Schroeder, K., Martsters, S.A., 
Koeppen, H., Ashkenazi. A ., and Kim, K.J. Isotype-dependent inhibition of 
tumor growth in vivo by monoclonal antibodies to death receptor 4. J. Immunol. 
166,4891-4898(2001). 

63. Pnllar.lf T F ^ F.rff. M.. and Ashkenazi. A . Direct stimulation of apoptotic 
signaling by soluble Apo2L/tumor necrosis factor-related apoptosis-inducing 
Ugand leads to selective killing of glioma cells. Clin. Cancer Res. 7, 1362-1369 
(2001). 

64. Wang, H., Marsters, S. A., Baker, T., Chan, B., Lee, W.P., Fu, L., Tumas, D., Yan, 
M., Dixit, V.M., * Ashkenazi. A ., and *Grewal, I.S. TACI-ligand interactions are 
required for T cell activation and collagen-induced arthritis in mice. Nature 
Immunol. 2, 632-631 (2001). 

65. Kischkei, F.C., Lawrence, D. A., Tinel, A., Virmani, A., Schow, P., Gazdar, A., 
Blenis, J., Amott, D., and Ashkenazi. A . Death receptor recruitment of 
endogenous caspase- 1 0 and apoptosis initiation in the absence of caspase-8. J. 
Biol. Chem. 216, 46639-46646 (2001). 

66. LeBlanc, H., Lawrence, D. A., Varfolomeev, E., Totpal, K., Morlan, J ., Schow, P., 
Fong, S., Schwall, R., Sinicropi, D., and Ashkenazi. A T umor cell resistance to 
death receptor induced apoptosis through mutational inactivation of the 
proapoptotitcBcl-2homologBax.A^afMreMerf. 8, 274-281 (2002). ■ 

67. Miller, K., Meng, G-, Liu, J., Hurst, A., Hsei, V., Wong, W-L., Ekert, R., 
Lawrence, D., Sherwood, S., DeForge, L., Gaudreault., Keller, G., Sliwkowski, 
M., Ashkenazi. A ., and Presta, L. Design, Construction, and analyses of 
multivalent antibodies. J. Immunol. 170, 4854-4861 (2003). 
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68. Varfolomeev, E., Kischkel, F., Martin, F., Wanh, H., Lawrence, D., Olsson, C, 
Tom, L., Erickson, S., French, D., Schow, P., Grewal, I. and Ashkenazi, A. 
Immune system development in APRIL knockout mice. Submitted. 

Review articles: 

1 . Ashkenazi. A.. Peralta, E., Winslow, J., Ramachandran, J., and Capon, D., J. 
Functional role of muscarinic acetylcholine receptor subtype diversity. Cold 
Spring Harbor Symposium on Quantitative Biology. LIII, 263-272 (1 988). 

2. Ashkenazi. A ., Peralta, E., Winslow, J., Ramachandran, J., and Capon, D. 
Functional diversity of muscarinic receptor subtypes in cellular signal 
transduction and growth. Trends Pharmacol. Sci. Dec Supplement, 12-21 (1989). 

3. Chamow, S., DuUege, A., Ammann, A., Kahn, J., Allen, D., Eichberg, J ., Byrn, 
R., Capon, D., Ward, R., and Ashkenazi. A . CD4 imlnunoadhesins in anti-HIV 

. therapy: new developments. J. Cancer Supplement 7, 69-72 (1992). 
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9. Chamow, S., and Ashkenazi. A . Immunoadhesins: Principles and Applications. 
Trends Biotech. U, 51-60 {\996): 

10. Ashkenazi. A ., and Chamow, S. M. Immunoadhesins as research tools and 
therapeutic agents. Curr. Opin. Immunol. 9, 195-200 (1997). 

11. Ashkenazi. A ., and Dixit, V. Death receptors: signaling and modulation. Science 
281, 1305-1308(1998). 

12. Ashkenazi. A ., and Dixit, V. Apoptosis control by death and decoy receptors. 
Curr.Opin. Cell. Biol. 11,255-260 (1999). 



9 



13. A^hkenazi. A . Chapters on Apo2I7TRA]L; DR4, DR5, DcRl, DcR2; and DcR3. 
Online Cytokine Handbook rwww.apnet.c oTn/cvtokinereference/). 

14. Ashkenazi, A . Targeting death and decoy receptors of the tumor necrosis factor 
superfamily. Nature Rev. Cancer 2, 420-430 (2002). 

15. LeBlanc, H. and Ashkenazi. A . Apoptosis signaling by Apo2I/rRAIL. Cell Death 
and Differentiation 10, 66^15 (2003). 

16. Ahnasan, A. and Ashkenazi. A . Apo2UTRAIL: apoptosis signaling, biology, and 
potential for cancer therapy. Cytokine and Growth Factor Reviews 14, 337-348 
(2003). 



Book: 

Antibody Fusion Proteins (Chamow, S., and Ashkenazi, A., eds., John Wiley and 
Sons Inc.) (1999). ... 



Talks: 

1 . Resistance of primary HIV isolates to CD4 is independent of CD4-gpl20 binding 
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Genome-wide Study of Gene Copy Numbers^ 
Transcripts, and Protein Levels in Pairs of 
Non-invasive and Invasive Human Transitional 
Cell Carcinomas* 

Torberi F. 0mtofti:§, Thomas ThykjaerU, Frec(eric M. Waldman||, Hans Wolf**, 
and Julio £• Cells** 



Gain and loss of chromosomal material is characteristic 
of bladder cancer, as well as malignant transformation in 
general. The consequences of these changes at both the 
transcription and translation levels Is at present unknown 
partly because of technical limitations. Here we have at- 
tempted to addre»ss this question in pairs* of non-invasive 
and Invasive human bladder tumors using a combination 
of technology that included comparative genomic hybrid- 



phenomenon at both the transcription and translation levels. 
High throughput anray studies of the breast cancer cell line 
BT474 has suggested that there Is a conelatlon between 
DNA copy numbers and gene expression in highly amplified 
areas (2), and studies of Individual genes In solid tumors 
have revealed a good correlation between gene dose and 
mRNA or protein levels in the case of c-ert>-62, cyciln dl, 
emsl, and N-myc (3-5). However, a high cycHn D1 protein 
expression has been observed without simultaneous am- 



ization, high density oligonucleotide array-based monitor- 
ing of transcript levels (5600 genes), and high resolution „^ ^ , , , ^ 
two-dimensional gel electrophoroslsfthe results showed^*^t«>n and a low level of c-myc copy number In- 
that there is a gene dosage effect 9»at ki some cases crease was obsen^ed without concomitant c-myc protein 
superimposes on other regulatory mechanisms. This ef- 
fect depended (p < 0«015) on the magnitude of the com- 
parative genomic hybridization change. In general (18 of 
23 cases], chromosomal areas with more than 2-fold gain 
of DNA showed a corresponding increase In mRNA tran- 
scripts. Areas with loss of DNA, on the other hand, 
showed either reduced or unaltered transcript levels;)Be- 
cause most proteins resolved by two-dimensional gels 
are unknown it was only possible to compare mRNA and 
protein alterations iixrelativefy few cases of well focused 
abundant proteins, ^ith few exceptions we found a good 
correlation (p < 0*00S) between transcript alterations and 
protein levels. The implications, as well as limitations, 
of the approach are discussed. Molecular A Cellular 
Proteomlcs 1:37-^45, SOOZ 



Aneuploidy Is a common feature of most human cancers 
(1), tHit little is known about the genome-wide effect of this 
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overexpression (6), 

in human bladder tumors, karyotyping, fluorescent In sftu 
hytMidization, and comparath^e genomic hybridization {CGHp 
have revealed chromosomal atienrations that seem to t>e 
characteristic of certain stages of disease progression, in the 
case of non-Invasive pTa transitionai cell carcinomas (TCCs), 
this includes loss of chromosome 9 or parts of K, as well as 
loss of Y In males. In minimally Invasive pTt TCCs, the fol- 
lowing alterations have been reported: 2q-, 11p-, 1q+, 
I1q13+, 17q+;and 20q+ (7-12). it has been suggested that 
these regions hart>or tumor suppressor genes and orico- 
genes; however, the large chromosomal areas Invoh/ed often 
contain many genes, making meartlngful predictions of the 
functional consequences of lasses and gains very difficult 

In this investigation we have combined genome-wide tech- 
nology for detecting genomic gains and losses (CGH) with 
gene expression profiling techniques (microanrays and pro- 
teomlcs) to determine the effect of gene copy number ori 
transcript and protein levels in pairs of non-lnvastve and In- 
vasive human k^ladder TGCs. 

EXPERIMENTAL PROCEDURES 

Materia/— Bladder tumor t>k>psies were sampled after Infonned 
cor>sent was obtained and after removal of tissue for routine pathol- 
ogy examination. By light microscopy tumor9.335 and 532 were 
staged by an experienced pathologist as pTa (superficial papillary), 

TTie abbreviations used are: CGH, conparative genomic hybrid- 
ization: TCC, transitional celt carcinoma: IjOH loss of heterozygosity; 
PA-FABP, psoriasis-associated fatty acid-bMing protein; 20, 
two-dimenstonaL 
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FiQ 1 DNA copy number and mRNA expression leveL Shown from teft to right are chromosome (ChrJ, CGH profiles, gene location and 
nmiJon level of specific genes, and overall expression level along the chromosome. A, expression of mRNA In Invasive tumor 733 as 
^^^^with the rwv4fivaslvo counterpart tumor 335. 6, expression of mRNA In invasive tumor 827 compared with the non-fcwasfve 
*^^rDart tumor 532 The average fluorescent signal ratio between tumor DNA and nonmal DNA Is shown alofig the length of the chromosome 
«m ti^^jW curve in the ratio profile represents a mean of four chromosomes and is swounded by thfn curves frxlicating one standard 
ntSfltten The centra/ verticai One prvken) indicates a ratio value of 1 (no change), and the vertical lines next to It (dotted) Indicate a ratio of 
oT^m md 2 0 iright m chromosomes where the non-lnvaslve tumor 335 used for comparison showed alterations In DNA content, the ratio 
nmfiteof that chrwnosome Is shown to the r^t of the Invash^e tumor profiler The cotored bars represents one gene each. Identified by the 
^'^^^r^mbers above the ban (the name of the gene can be seen at www,MDLDK/sdata.htmO. The bars indicate the purported location of 
II«nSrand the cotors hdicate the expression level of the gene m the Invasive tumor compared with the noa-4nvash« counterpart; >2-fold 
^!nr^a totocW >2-fold decrease (blue), no signlficanl change (orange). The bar to the fer nSff/if. entitled &pre$slon shows the resulting change 
r^®^!^^ the chromosome: the cotors Indicate that at least half of the genes were up-regulalad (Wack), at least half of the genes 
H^iyy!^^lated (Wue) or more than half of the genes are unchanged (omnige). If a gwie was absent hi one of the samples and present In 
regarded as more than a 2-fold change. A 2-fold level was chosen as this corresponded to one standaird deviation In a double 
^ennlrurtlon of --1800 geries. Centromeres and heterochforn^ 

grade I and respective, tumors 733 and 827 were staged as pTI 
Onvaslve Into submucosa). 733 was staged as solid, and 827 was 
staaed as papillary. t)oth grade III. 

mflWA Pr^tfon-Tlssue biopsies, obtained fresh from surgery, 
were embedded immediately In a sodium-guanldinium thiocyanate 
solution and stored at -80 'C. Total RNA was isolated usihg the 
B RNA isolation method (WAK-Chemie Medical GMBH). 
poly(A)* RNA was Isolated by an ollgo(dT) selection step (Ollgotex 

cflWA*Pre/^^tort^^ m9 °' mRNA was used as startir)g material. 
The fir^ and second strand cDNA synthesis was performed using the 
SuperScripKB* choice system (Invitrogen) according to the manufac- 
turer's instructions but using an ofigo{dD primer containing a T7 RfMA 
DOlymerase binding site. Latwied cRNA was prepared using the ME- 
GAscrp® h vitro transcription Wt (AmWon); Biotfn-labeled CTP and 



UTP (Enzo) was used, together with unlabeled NTPs in the reaction. 
Foltowlr^g the in vitro transcription reacUoa the unincorporated nu* 
deotldes were renraved using RNeasy columns (Qiagen). 

Array Hybridization and Scanning— Arra:^ hybriklization and scan- 
ning was modified from a previous method (13). 10 |ig of cRNA was 
fragmented at 94 MS for 35 min in t>uffer containing 40 mM Tris 
acetate, pH 8.1 , 1 00 mM KOAc. 30 mM MgOAc. Prior to hybridization, 
the fragmented cRNA in a 6X SSPE-T hybridizatior) buffer (1 m NaCl. 
10 mM Trfs, pH 7^, 0.00596 Tritori), was heated to 95 *C for 5 min, 
subsequently cooled to 40 *C, and loaded onto the Affymetrix probe 
array cartridge. The prot>e anray was then lncut)ated for 16 h at 40 *C 
at constant rotation (60 rpm). The probe array was exposed to 10 
washes in 6x SSPE-T at 25 *C followed by 4 washes In 0.5x SSPE-T 
at 50 *C. The bioUnylated cRNA was stained with a streptavidin- 
phycoerythrin conjugate, 10 /tg/ml (Molecular Probe^ in 6x SSPE-T 
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for 30 min at 25 -C followdd by 1 0 washes In 6x SSPE-T at 25 'C, The 
probe arrays were scanned at 560 nm using a confocal laser scanning 
microscope (made lor AffymetrU by Hewtett-Packard). The readings 
from the quantitative scanning were anaJyzed by Affymetrix gene 
expression analysis software, 

Mhuossteme AnatyslS'-f^^crxssaXme Analysis was perfwrned as 
described previously (14). MIcrosatelliles were selected by use of 
wwwjwbljiim.nlh.gov/genemap98, and primer sequences were ob- 
tained from the genome data base al www.gdb.org. DNA was extracted 
from tumor and Wood and ampfifled by PGR In a volume of 20 Ml for 35 
cycles Ihe ampllcons weredenatured and electrophoresed for 3 h In an 
ABI Prtem 377. Data were coHecled in the Gene Scan program for 
fragment analysis. Loss of heterozygosity was defined as less than 33% 
of one aneie delected in tunor amplicons compared with Wood. 

Pfoteomfc Anaiysis-tCCs were mipced Into smalj pieces and 
homogenized in a small glass homogenlzer in 0.5 ml of lysis solution. 
Samples were stored at -20 unttt use. The procedure for 2D get 
electrophoresis has been descrft)ed in detail elsewhere (15. IQ. Gels 
wera stained with silver nitrate and/or Coomassle Brilliant Blue. Pro- 
teins were Weniffied by a combination of procedures that Included 
microsequendng. mass spectrometry, twoHjimensional gel Western 
immunoWotting, and comparison with the master two-dimensional gel 
image of human keratinoc/te proteins; see biobase.dk/<^i-birVcefis. 

CGH* Hybridization of differentially lat>eled tumor and normal DIMA 
to normal metaphase chromosomes was performed as descrtk>ed 
previously (1(3). Ruorescein-labeled tumor Df^ (200 ng). Texas Red- 



labeled reference OHA (200 ng), arxj human Cot-1 DNA (20 fig) were 
denatured at 37 *C for 5 nnln and applied to denatured r>onnal met- 
aphase sfides. Hyt>ridizat3on was at 37 *C for 2 days. After wasMng, 
the slides were counterstained wf^ 0.15 ftg/ml 4,6-diamldno-2-phe* 
nylindole in an anti-fade sotutlon. A second hybridization was per- 
fonmed for all tumor samples using fiuoresceiiHabeled refersnce DNA 
arxi Texas Red-labeled tumor DNA (inverse labeling} to confirm the 
aberrations detected during the initial hybrtcflzatfon. Each CGH ex- 
periment also included a normal control hybridizalion usirtg fluores- 
cein- and Texas Red-labeled rtormal Df^lA. Digital image analysis was 
used to identify chromosomal regions with abnormal fluorescence 
ratios, indicating regions of DNA gains and losses. The average 
green:red fluorescence intensity ratio profSes were calculated using 
four images of each chromosome (eight chromoson>es totaO With 
normalization of the greert:red fluorescence Intensity ratio for the 
entire metaphase and background correction. Chromosome Identifi- 
cation was performed t>ased on 4,6-dtamidino-2-phenytirKtole tsand- 
ing patterns. Only images showing uniform high intensity fluores- 
cence with minimal backgrourKi staining were ana^fzed. AU 
centrom^es, p arms of acrocentric chromosomes, and heterochro- 
matic regions were excluded from the artalysis. 

RESULTS 

Comparath/e Genomic Hyi^dizatSon-^lhe CGH analysis 
identified a number of chromosomal gains and losses in the 



Molecular & Cellular Proteomics f « 1 39 



(■ 



Gene Copy Numbers, Transcripts, and Protein Levels 



Table I 



Correladon between an&Btions detected by (XSH and 

Top, CGH used as Independent variable (if CGH alteration - what expression ratio was found); bottofn, aHered express used i 
independent variable (ti expresston aHeration - what CGH deviation was found). 



CGH alterations 



Tumor 733 vs. 335 
Expression ctenge dusters 



Concordance CGH alterations 



Tumor 827^ vs, 532 
Expression change clusters 



Concordance 



13 Gain 



10 Loss 



10 Up-fegutation 

0 Down-ragutalion 

3 No change i' 

1 Up^regulatioh 

5 Down-regulation 

4 No change 



77% 



5096 



10 Gain 8 Up-rsgUatlon 

0 OownrTsgulation 
2 No change 

12 Loss 3 Up-rsgulation 

2 Down regutartion 
7 r^o change 



Expression change clusters 



TunK>r733vs.335 
CGH alterations 



Concordance Exprsssion chartge clusters 



Tumor a27 vs. 532 
CGH alterations 



80% 
17% 

Concordance 



16 Up-reguiation 
21 Down-regulation 
15 No change 



llOatn 

2 Loss 

3 No change 
1 Gain 

8 Loss 

12No change 
3 Gain 
3 Loss 

9 No change 



69% 
38% 
60% 



17 U^MegutatkMi 
9 Down-regulation 
21 No change 



10 Gain 

5 Loss 

2 No ch&nge 
OGain 

3 Loss 

6 No chainge 
1 Gain 

3 Loss 

17Noc^^ange 



59% 
33% 
81% 



two invasive tumors {stage pT1. TCCs 733 and 827), whereas 
the two non-invasive papillomas (stage pTa, TCX^s 335 and 
532) showed only 9p-. 9q22-q33-, and and 7+. 9q-, 
and respectively. Both invasive tumors showed chaniges 
(1q22-24+, 2ql4.1-qter-, 3q12-q13^-. 6q12-q22-, 
9q34+, 11q12-q13+,M7+, and 20q11,2-q12+) that are typ- 
ical for their disease stage, as well as additional alterations, 
some of which are shown In Fig, 1. Areas with gains and 
loi^es deviated from the normal copy number to some extent, 
and the average numerical deviation from rxmnal was 0.4-fotd 
In the case of TCC 733 and 0.3-fold for TCC 827. The largest 
changes, amounting to at least a doubling of chromosomal 
content, were observed at 1q23 In TCC 733 (Rg. ^A) and 
20q12 In TCC 827 (Rg. IS). 

mRNA Expri^lon in Reiation to DMA Copy Numbet'-The 
mRfviA levels from the two Invasive tumors (TCCs 827 and 
733) were compared with the two non-Invasive counterparts 
(TCCs 532 and 335). This was done In two separate experi- 
ments in which we compared TCCs 733 to 335 and 827 to 
532, respecth^ety. using two different scaling settings for the 
anays to rule out scaling as a confounding parameter. Ap- 
proximately 1.800 genes that yielded a signal on the arrays 
were searched in the Unigene and Genemap data bases for 
chromosomal location, and those with a known location 
(1096) were plotted as bars covering their purported locus. In 
that way it possible to construct a graphic presentation of 
DNA copy number and relative mRNA levels along the Indi- 
vidual chroniosomes (Rg. 1). 

For each mRNA a ratio was calculated between the level in 
the invasive versus the non-invasive counterpart Bars, which 
represent chromosomal location of a gene, were color-coded 
according to the expresston ratio, and only differences larger 



than 2-fold were regarded as Irrformath/e (Rg, 1). The density 
of genes along the chromosomes varied, and areas contain- 
ing only one gene were excluded frorh the calculations. The 
resolution of the CGH method Is very low, and some of the 
outlier data may be because of the fact that the botmdaries of 
the chromosoHDal aberrations are not kntiym at high resolutioa 
Two sets of calculations were nnade from the data. For the 
first set we used CGH alterations as the Independent variable 
and estimated the frequency of compression alterations In these 
chromosomal areas. In general, areas with a strong gain of 
chromosomal material contained a cluster of genes havirfg 
Increased mRNA expression. For example, both chrorno- 
somes 1q21-q25, 2p and 9q, showed a relative gain of more 
than 100% in DNA copy number that was accompanied by 
increased mRNA expression levels in the two tumor pairs (Rg. 
1). in most cases, chromosomal gains detected by CGH were 
accompanied by an increased level of transcripts in both 
TCCs 733 (77%) and 827 (80%) (Table I, top). Chromosomal 
losses, on the other hand, were riot accompanied by de- 
creased expression In several caees, and were often regis- 
tered as having unaltered RNA levels (Table I, fop). The Inabil- 
ity to detect RNA expression changes In these cases was not 
because of fewer genes mapping to the lost regions (data not 
shown). 

In the second set of calculations we selected expression 
alterations above 2-fold as the independent variable and es- 
timated the frequency of CGH alterations in these areas. As 
above, we found that Increased transcript expression corre- 
lated with gain of chromosomal material (TCC 733. 69% and 
TCC 827, 59%), vrhereas reduced expression was often de- 
tected in areas with unaltered CGH ratios (Table I. bottom). 
Furthermore, as a control we looked at areas with no atter- 
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Fio. 2. Correlation between maximum CGH aberration and the ability to detect expression change by oligonucleotide anr^ 
monitoring. The aberration Is shown as a numerical -fold change in ratio between invasive tun>ors 827 (A) and 733 (4 ) and their norv4nvasive 
counterparts 532 and 335. The expression change was taken from the Expression line to the r^t In Fig. 1. which depicts the resulting 
expression charge for a given chromosomal region. At least half of the mRNAs from a given region have to be either up- or down-regulated 
to be scored as an expression change. All chromosomal arms in which the CGH ratio plus or minus one standard de\^ion was outside the 
ratio value of one were included. 



ation In expression. No alteration was detected by CGH In 
most of these areas (TCC 733, 60% and TCC 827, 81 %; see 
Table I, bottom)* Beceuse the at>iltty to observe reduced or 
Increased mRNA expression clustering to a certain chromo- 
somal area clearly reflected the extent of copy number 
changes, we plotted the maximum CGH aberrations In the 
regior^s showing CGH changes against the ability to detect a 
change in mRMA expression as monitored by the oligonucleo- 
tide arrays (Rg. 2)(fi)r both tumors TCC 733 ft? < 0.01 §) and 
TCC 827 (p < 0.00003) a highly significant conrelation was 
observed between the level of CGH ratio change j(reflecting 
the DNA copy number) and aiterations detected by the array 
based technology (Fig. ^ Similar data were obtained when 
areas with altered expression were used as independent vari- 
ables. These areas con^elated best with CGH when the CGH 
ratio deviated 1.6- to 2.0*fold (Table I, bottom) but mostly did 
not at lower CGH deviations. These data probably reflect that 
loss of an allele may only lead to a 50% reduction In expres- 
sion level, which is at the cut-off point for detection of expres- 
sion alterations. Gain of chromosomal material can occur to a 
much larger extent 

Microsatelltte-based Detecpon of Minor Areas o/ Loss- 
es— In TCC 733, several chromosomal areas exhibiting DMA 
amplification were preceded or followed by areas with a nor- 
mal CGH but reduced mRNA expression (see Rg. 1, TCC 733 
chromosome 1q32, 2p21, and 7q21 and q32, 9q34, and 
10q22}. To determine whether these results were because of 
undetected loss of chromosomal material In these regions or 



t>ecause of other non-structural mechanlsnDS regulating tran« 
scriptlon, we examined two microsatellites positioried at chro- 
mosome 1q25-32 and two at chromosome 2p22. Loss of 
heterozygosity (LOH) was found at both 1q25 and at 2p22 
indicating that minor deleted areas were not detected with the 
resolution of CGH (Rg. 3). Addttionalty, chromosome 2p In 
TCC 733 showed a CGH pattern of gain/ho change/gain of 
DNA that conrelated with transcript Increase/decrease/in- 
crease. Thus, for the areas showing increased expresdon 
there w^ a correlation with the DNA copy number afterations 
(Rg. 1^). As Indicated above, the mRNA decrease ol>served in 
the middle of the chromosomal gain was because of LOH, 
implying that one of the mechanisms for mRNA down-regu- 
lation may be regions that have undergone smaller losses of 
chromosomal material. However, tills cannot be detected with 
the resolution of the CGH method. 

In both TCC 733 and TCC 827, the telomeric end of chro- 
mosome 11 p showed a nomtal ratio In the CGH analysis; 
however, clusters of five and three genes, respecUvely, lost 
their expression. Two microsatellites (D11S1760, D11S922) 
positioned close to MUC2. IGF2, and cathepsln 0 indicated 
LOH as the most llicely mechanism t^ehlnd the loss of expres- 
sion (data not shovym). 

A reduced expression of mRNA observed in TCC 733 at 
chromosomes dq24, 11p11, 12p12.2. 12q21.1, and 16q24 
and in TCC 827 at chromosome 11p15.5, 12p11, 16q11.2, 
and 18q12 was also examined for chrpmosonnal losses using 
microsatellites positioned as close as possible to the gene loci 
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Fio. 3. Microsateltite analysis of loss of heterozygosKy. Tumor 
733 showing loss of heterozygosity at chromosome 1q25. detected 
(a) by D1S215 dose to Hu dass I Nstocompatlblllty antigen (gene 
number 38 in Fig. 1), (b) by D1S2735 close to cathepsin E (gene 
number 41 in Rg. 1). and (c) at chromosome 2p23 by 02S2251 dose 
to general 0-spectrin (gene number 1 1 on Fig. 1) and o( (a) tumor 827 
showing loss of heterozygosfty at chromosonne 18q12 by S18S1118 
dose to mitochondrial 3-oxoacyf-coenzyme A thiolase (gene numt>er 
12 In Fig. 1). The upper curves show the eleclropherogram obtained 
from normal ONA from leukocytes (AO, and the /ower ctirves show the 
electropherograiT) from tumor DMA (7). In all cases one aDele is 
partfetfly tost In the tumor amplioon- 

showing reduced mRN A transcripts. Only the microsatellite 
positioned at 18q12 showed LOH (Rg. 3], suggesting that 
transcriptional down-regulation of genes In the other regions 
may be controlled by other mechanisms. 

Relatfon between Changes in mRNA and Protein Levefe— 
2D-PAGE analysis, in combination with Coomassie Brilliant 
Blue and/or silver staining, was earned out on all four tumors 
using fresh biopsy material, 40 well resolved abundant known 
proteins migrating In areas away from the edges of the pH 
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FiQ. 4. Correlation between protein levels at Judged by 20- 
PAGE and transcript ratio. For comparison proteins were divkled In 
three groups, unaltered In level or up- or down-regulated (Pxvtentaf 
. axis). The mRNA ratio as detenmined tay digonudeotkie arrays was 
plotted for each gene {verdcai axis), mRNAs that were scored as 
present In k>oth tumors used for tt>e ratio calculation; ^ mRNAs that 
were scored as absent In the Invasive tumors (atong hoftiontataxls)or 
as absent in norHrivasive reference (top of figure). Two different 
scalings were used to exdude scaling as a confounder, TCCs 827 
arKl 532 (A^ were scaled wfth iMck^cund suppression, aiKJ TCCs 
733 and 335 (•O) were scaled without sippresslort Both comparf- 
sons showed highly significant (p < 0.005) differences in mRNA ratios 
t>etween the groups. Proteins shown were as foUov^ Group A (from 
iefti, phosphogtucomutase 1 , glutathtor»e transferase class |i rHimt>er 
4, fatty add-t>inding protein homdogue, cytokeratin 15, and cyto- 
keratin 13; e (from teff), fatty add-bindlrig protein homdogue, 28-kOa 
heat shocl< prbtda cytdceratin 13, and calcycUn; C<from left^, a-eno- 
lase, hnRNP 81, 28-kDa heat shock protein, 14-0-3-^ and 
pre-mRNA spUdng factor; 0, mesothelial keratin K7 (type 11); £ (from 
top), glutathk>f>e S-transferase-ir and mesothefial keratin K7 (type II); 
F (from top and left), adenyiyi cydase-assodated protein. E-cadherin, 
keratin 19, calgizzarin, phosphoglycerate mutase, annexin (V. cy- 
toskeletal r^ctin, hnRNP AI, integral membrane protein calnexin 
(IP90), hnRNP H, k>rain-type dathrin tight chain-a, hnRNP F, 70-fcOa 
heat shock protein, heterogeneous nudear ribonudeoprotdn A/B. 
tran$latk>nalV contrdied tumor protein, liver glyceraldehyde-3-phos- 
pfiate dehydrogenase, keratin 8, aklehyde reductase, and Na,K- 
ATPase 0-1 8ut>untt; G, (from top and ieffj, TCP20, calgizzarin, 70- 
kDa heat shock protein, calnexin, hnRNP H, cytokeratirt IS, ATP 
synthase, keratin 1 9, trk>sephosphate isomerase, hnRNP liver glyc- 
eraldehyde-3-phosphata3e dehydrogenase, glutathk>ne S-transfer- 
ase-w, and keratin 8; H (from left), plasma gdsoOn, autoant^ ca>- 
raticunn, thioredoxin, and NAD+ -dependent 15 hydroxyprbstagtandn 
dehydrogenase; f (from top), prolyl 4-r^roxylase /^-sutxmit, cyto- 
keratin 20, cytokeratin 17, prohibftlon, and fnictose 1,6-biphos- 
phatase; J annaxin »; K, annoxfri IV; L (from top and tefl), 90-kDa heat 
shock protein, prdyl 4-hydrQxylase 0-sutxinit, a-endase, GRP 78, 
cyctophilia and cofdin. 

gradient, and having a known chromosomal bcation, were 
selected for analysis in the TCC pair 827/532. Proteins were 
identified by a comt^natlon of methods (see ''Experimental 
Procedures"). In general there was a high^ significant corre- 
lation (p < 0.005) between mRNA and protein alterations (Rg. 
4). Only one gene showed disagreement between transcript 
alteration and protein alteratton. Except for a group of cyto- 
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Fig. 5. Comparison of protein and transcript lavets In invasive 
and non-invasive TCCs. The upper part of the figure shows a 2D gel 
ifeft) and the digonudeotide array (h^O ofTCC Saa-Theretfrecfian- 
gtes on the upper gel highlight the areas that are compared below. 
Identical areas of 20 gels of TCX^s 532 and 827 are shown below. 
Clearly, cytokeratins 13 and 15 are strongly dowrweguiated k\ TCC 
827 (fed annotatlcn]. The tile on the array contalr>ing prot>es for 
cytokeratin 15 is enlarged below the array (rec/ arrow) from TCC 532 
and is compared with TCC 827, The upper row of squares (n each tile 
corresponds to perfect match probes; the hwer row corresppr>ds to 
mismatch prot>es containing a mutation (used for correction for un- 
speciflc binding]. Absence of signal Is depicted as black, and the 
higher the signal the fighter the color. A high transcript level was 
detected in TCC 532 (6151 units) whereas a much lower level was 
detected in TCC 827 (absence of signals). For cytokeratin 13, a high 
transcript level was also present in TCC 53i2 (15659 units), and a 
much k>wer level was present in TCC 827 (623 units). The 2D gets at 
the bottom of the figure (te/© show levels of PA-FABP and adipocyte- 
FABP In TCCs 335 and 733 (invasive), respectively. Both proteins are 
down-regulated in the invasive tumor. To the right we show the array 
tiles for the PA-FAB P transcript A medium transcript level was de- 
tected In the case of TCC 335 (1277 units) whereas veiy low levels 
were detected in TCC 733 (166 units). /£F, isoelectric focusing. 



keratins encoded by genes on chr-onrK>some 17 (Rg, 5) the 
analyzed proteins did not belong to a particular fami 26 well 
focused proteins whose genes had a know d)romosomal 
location were detected In TCCs 733 and 335, and of these 19 
conrelated (p < 0.005) with the mRM A changes detected using 
the anays (Rg. 4). For example, PA-FABP was higWy ex- 
pressed in the non-invasive TCC 335 but k>st In the Invash/e 
counterpart (TCC 733; see Rg. 5>. The smafler numt>er of 
proteins detected In both 733 and 336 was because of the 
smaller size of the biopsies that weie avall£d3le. 

1 1 Chromosomal regions where CGH showed aberrations 
that conresponded to the changes In trartscrlpt levels also 
showed corresponding changes In the protein level (Table IQ. 
These regions Included genes that encode proteins that are 
found to be frequently altered in t>ladder carreer, nam^ 
cytokeratins 17 and 20, anhexins li and IV, and the fatty 
acid-binding proteins PA-FABP and FBP1, Four of these pro- 
teins were encoded by genes in chromosome 17q. a fre- 
quently amplified chromosomal area In Invasive bladder 
cancers* 

DISCUSSION. 

f^ost human cancers have abriormal DMA content, having 
lost some chromosomal parts and gained others. The present 
study provides some evidence as to the effect of these gains 
and losses on gene exp>resslon in two pairs of norHnvasive 
and invasive TCCs lising high throughput expresston arrays 
and proteomtos, In combination vrlth CGH. In general, the 
results showed that there Is a clear IrKilvkjual regulation of the 
mRNA expression of single genes, which In some cases was 
superimposed by a DNA copy number effect In most cases, 
genes located in chromosomal areas with gains often exhlt>- 
ited increased mRNA ex|3ression, whereas areas showing 
losses showed either no change or a reduced mRNA expres- 
sion. The latter might be because of the fact that k>sses most 
often are restricted to loss of one allele. arKi the cut-ofT point 
for detection of expression alterations was a 2-fold change, 
thus being at the t>order of detection. In several cases, how- 
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^ In cases where the corresponding alterations wera found in t^oth TCCs 827 and 733 these are shown as 827/733. 
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ever, an increase or decrease in DNA copy number was 
associated with de novo occurrence or complete loss of tran- 
script* respectively. Some of these transcripts could not be 
detected in the non-invash^e tumor but were present at rela- 
tively high levels In areas with DMA amplifications in the inva- 
sive tumors (ag. In TCC 733 transcript from cellular llgand of 
annexin II gene (chromosome 1q21) from absent to 2670 
arbitrary units; In TCC 827 transcript from small proiine-rich 
protein. 1 gene (chromosome 1ql2-<|21.1) from absent to 
1326 arbitrary units). It may be anticipated fmm these data 
that significant clustering of genes with an Increased expres- 
sion to a certain chromosdnDal area indicates an increased 
likelihood of gain of chromosome material In this area." 

Considering the many possible regulatory mechanisms act- 
ing at the level of transcription, it seems striking that the gene 
dose effects were so dearly detectable in gained areas. One 
hypottietk:al explanation may lie In the loss of controlled 
methylatlpn In tumor cells (17-1$). Thus, It niay be possible 
that In chromosomes with Increased DMA copy numbers two 
or more alleles could be demethyiated simultaneously leading 
to a Ngher transcription level, whereas In chronru>5omes with 
losses the remaining allele could be partly nr>ethylated, timing 
off the process (20, 21). A recent report has documented a 
ploldy regulation of gene expression in yeast, but In this case all 
the genes were prasemt in the same ratk) (2^, a situation that Is 
not analogous to that of canc^ cells, which show marked 
chromosomal abenations, as well as gene, dosage effects. 

Several CGH studies , of bladder cancer have shown that 
some chromosomal aberrations are cornmon at certain 
stages of disease progression, often occurring in more than 1 
of 3 tumors. In pTa tumors, these include 9p-, 9q-, 1 q +. Y- 
(2, 6), and In pT1 tunriors, 2q-,11p", 11q-, 1q+, 5p+, 8q+. 
17q+, and 20q+ (2-4, 6, 7). The pTa tumors studied here 
showed similar ab>errations such as 9p- and 9q22-q33- and 
9q^ and respectively. Likewise, the two minimal invasive 
pTI tumors showed aberrations that are commonly seen at 
that stage, and TCC 827 had a remarkable resemblance to the 
commonly seen pattern of losses and gains, such as 1q22-24 
amplification (seen in both tumors), 1 1q14-q22 loss, the latter 
often linked to 17 q+ (both tumors), and 1q+ and 9p-, often 
linked to 20q+ and 11 q13+ (both tumors) (7-9). These ob- 
servations indicate that the pairs of tunruxs used In this study 
exhibit chromosomal changes observed in many tumors, and 
therefore the findings could be of general Importarice for 
bladder cancer. 

Considering that the nnapping resolution of CGH is of about 
20 megabases it Is only possible to get a crude picture of 
chromosomal instability using this technique. Occasionally, 
we observed reduced transcript levels dose to or inside re- 
gions with Increased copy numt^ers. Analysis of these regions 
by positioning heterozygous microsateWtes as close as pos- 
sible to the locus showing reduced gene expression revealed 
loss of heterozygosity in several cases. It seems likety that 
multiple and different events occur along each chromosomal 



arm and that the use of cDNA microarrays for analysis of DNA 
copy number changes will reach a resolution that can resoh^e 
these changes, as has recently been proposed (2). The outlier 
data were not more frequent at the boundaries of the CGH 
abenations. At present we do hot icnow the mechanism be- 
hind chromosomal aneuptoldy and cannot predkH: whether 
chromosomal gains will be transcribed to a larger extent than 
the two native alleles. A mechanism as genetic Imprintir^ has 
an impact on the expression level In nprrhaJ celts and is often 
reduced In tumors. Howev^, the relation between Imprinting 
and gain of cfvomosomal material is not known. 

We regard It as a strengtti of this irivestigation that we were 
able to compare invasive tumors to benign turners rather than 
to normal urothelium, as the tumors studied were blolo^cally 
very close and probably may represent successive steps in 
the progression of bladder cancer. Despite the limited amount 
of fresh tissue available it was possible to apply three different 
state of the art methods. The observed conrelatk>n between 
DNA copy number and mRNA expression is remarkable when 
or>e considers that different pieces of the tumor bk}psies were 
used for the cpferent sets of experiments. This Indicate that 
bladder tumors are relativety homogenous, a notk>n recently 
supported by CGH and LOH data lhaX showed a remarkable 
similarity even between tumors and distant metastasis (10, 23). 

In the few cases analyzed, mRNA and protein levels 
showed a striking correspondence although In some cases 
we found discrepancies that may be attributed to translational 
regulation, post-translatlonal processing, protein degrada- 
tion, or a combination of these. Some transcripts belor^ to 
undertranslated mRNA pools, which are associated with few 
translationally inactive ribosonles; these pods, however, 
seem to bB rare (24). Protein degradation, for example, may 
be very Important in thp case of polypeptides with a short 
haJf-life (e.^. signaling proteins). A poor conrelation between 
mRNA and protein levels was found in Ih^ cells as deter- 
mined by arrays and 2D-PAGE (ISj, and a nrxxierate conela- 
tloQ was recently reported by Ideker et af. (26) In yeast 
(Tnterestingly, our study revealed a nuich better correlation 
between gained chromosomal areas and Increased mRNA 
levels than between loss of chronnosomal areas and reduced 
mRNA levels. In general, the level of CGH change determined 
the ability to detect a change in trariscripC) One possible 
explanation could be that by losing one allele the chainge In 
mRNA level Is not so dramatic as compared with gain of 
material, which can be rather unlimited and may lead to a 
severatfold increase In gene copy numtm resulting in a much 
higher Impact on transcript level. The latter would fc>e much 
easier to detect on the expression arrays as the cut-off point 
was placed at a 2-fold level so as not to be biased by noise on 
the array. Construction of arrays with a better signal to noise 
ratio may In the future altow detection of lesser than 2-fold 
alterations in transcript levels, a feature that may facilitate the 
analysis of the effect of bss of chromosomal areas on tran- 
script levels. 
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In eleven cases we found a significant correlatioh between 
DMA copy number, mRNA expression, and protein level. Four 
of these proteins were encoded by genes located at a fre- 
quently amplified area in chromosome 17q. Whether DMA 
copy numt>er is one of the mechanisms behind alteration of 
these eleven proteins Is at present unknown and will have to 
be proved by other methods using a larger number of sam- 
ples. One factor making such studies complicated Is the targe 
extent of protein modification that occurs after translation, 
requiring Immunoidentiflcation and/or mass spectrometry to 
con^ectly identify the proteins in the gels. 

In conclusion, the results presented In this study exemplify 
the large body of knowledge that may be possible to gather In 
the future t>y combining state of the art techniques that follow 
the pathway from DMA to protein (26). Here, we used a tradi- 
tional chromosomal CGH method, but in the future high reso- 
lution CGH based on mk^roarrays with many thousand radiation 
hybrid-mapped genes v»/lll Increase the resolution and infbrma- 
tlon derived from these types of experiments p). Combined with 
expressk>n arrays analyzlrig transcripts derived from genes with 
known locations, and 2D gel analysis to obtain Information at 
the post-translational level, a clearer and more developed un- 
derstanding of the tumor genome win be forthcoming. 
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ABSTRACT 

Genetic changes underlie tumor progression and may lead to cancer- 
specific expression of critical genes. Over 1100 pubUcations have de- 
scribed the use of comparative genomic hybridization (CGH) to analyze 
the pattern of copy number alterations in cancer, but very few of the genes 
affected are known. Here, we performed high-resolution CGH analysis on 
cDNA microarrays in breast cancer and directly compared copy number ^ 
and mRNA eipression leveU of 13,824 genes to quantitatc the Impact of 
genomic changes on gene e^qpression. We identified and mapped the 
boundaries of 24 independent amplicons, ranging in sixe from 0*2 to 32 
Mb. Throughout the genome, both high- and low-level copy number 
changes had a substantial impact on gene expression, with 44% of the 
highly ampUfied genes showing overexprcssion and 10.5% of the highly 
overcxpressed genes being aropUfied. Statistical analysis with random 
permutation tests identified 270 genes whose expression levels across 14 
samples were systematically attribuUble to gene ampUfication. These 
included roost previously described amplified genes in breast cancer and 
many novel targets for genomic alterations, including the H0XB7 gene, 
the presence of which in a novel amplicon at 17q213 was validated In 
10.2% of primary breast cancers and associated with poor patient prog- 
nosis, in conclusion, CGH on cDNA microarrays revealed hundreds of 
novel genes whose overexpression Is attributable to gene amplification. 
These genes may provide insights to the clonal evolution and progression 
of breast cancer and highlight promising therapeutic targets. 

INTRODUCTION 

Gene expression patterns revealed by cDNA niicroaTrays have 
facilitated classification of cancers into biologically distinct catego- 
ries, some of which may explain the clinical behavior of the tumors 
(1-6). Despite this progress in diagnostic classification, the molecular 
mechanisms underlying gene expression patterns in cancer have re- 
mained elusive, and the utility of gene expression profiling in the 
identification of specific therapeutic targets remains limited. 

Accumulation of genetic defects is thought to underiie the clonal 
evolution of cancer. Identification of the genes that mediate the effects 
of genetic changes may be important by highlighting transcripts that 
are actively involved in tumor progression. Such transcripts and their 
encoded proteins would be ideal targets for anticancer therapies, as 
demonstrated by the clinical success of new therapies against ampli- 
fied oncogenes, such as ERBB2 and EGFR (7. 8), in breast cancer and 
other solid tumors. Besides amplifications of kiwwn oncogenes, over 
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Copy number ratio 




Expression ratio 

Fig. 1. Intact of gene copy number on global gene expression Icveb. A. percentage of 
over- and undtraprcsscd genes {X oxis) according to copy number ratios {X axis) 
Threshold values used for over- and undcrexpression were >2.I84 (global upper 7% or 
the cDN A ratios) and <0.4826 (global lower 7% of the expression ratios). B. percentage 
of amplified and deleted gcocs according to cxpicssion rarios. Threshold values for 
araplificatwn and deletion were > U and <0.7. 



20 recurrent regions of DNA amplification have been mapped in 
breast cancer by CGH^ (9. 10). However, these amplicons are often 
large and poorly defined, and their impact on gene expression remains 
unknown. 

We hypothesized that genome-wide identification of those gene 
expression changes that are attributable to underlying gene copy 
number alterations would highlight transcripts that arc actively in- 
volved in the causation Or maintenance of the malignant phenotype. 
To identify such transcripts, we applied a combinaUon of cDNA and 
CGH microarrays to: {a) determine the global impact that gene copy 
number variation plays in breast cancer development and progression; 
and (6) identify and characterize those genes whose mRNA cxpres- 



^ The abbreviations used are: CGH, comparative genomic hybridization; RSH. fluo- 
rescence Ml situ hybridization; RT-PCR, reverse transcription-PCR. 
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sion is most significantly associated with amplification of the corre- 
sponding genomic tempiaie. 

MATERIALS AND METHODS 

Breast Cancer CeJI Uries. Fourteen breast cancer ceU lines (BT-20, BT- 
474, HCC1428, Hs578t. MCF7, MDA-361, MDA-436, MDA-453, MDA^68, 
SKBR-3. T-47D, UACCS12. 2R-75-U and ZR*75-3Q) were obtained from the 
American Type Culnirc Collection (Manassas, VA). Cells were gro^n under 
recommended culture conditions. Genomic DNA and niRNA were isolated 
using standard protocols. 

Copy Number and Expression Analyses by cDNA Microarrays, The 
preparation and printing of die 13.824 cDNA clones on glass sHdcs were 
performed as described (1 1-13). Of these clones. 244 represented uncharac- 
terized «q>ressed sequence togs, and the rcmamder corresponded to known 
genes. CGH experiments on cDNA microanays were done as described (14, 
15). Briefly, 20 jxg of genomic DNA from breast cancer cell lines and normal 
human WBCs were digested for 14-18 h with Ahl and Rsal (Life Technol- 
ogies, Inc.. Rockville. MD) and purified by phenol/chloroform extraction. Six 
of digested cell line DNAs were labeled with Cy3-dUTP (Amcrsluim 
Pharmacia) and normal DNA with Cy5-dUTP (Amershara Pharmacia) using 
the Bioprime Labeling kit (Life Technologies, Inc.). Hybridization (14. 15) and 
posthybridizalion washes (13) were done as described. For the expression 
analyses, a standard reference (Universal Human Reference RNA; Siratagcne, 
U Jolla'CA) was used in all experiments. Forty /ig of reference RNA were 
labeled with Cy3-dUTP and 3.5 fig of test mRNA with Cy5-dUTP, and the 
labeled cDNAs were hybridized on microarrays as described (1 3, 1 5). For both 
microarray analyses, a laser confocal scanner (Agilent Technologies, Palo 
Alto. CA) was used to measure the fluorescence intcnsiUcs at the target 
locations using the DE ARRAY software (16). AAer background subtraction, 
average intensities at each clone in the lest hybridization were divided by the 
average intensity of the corresponding done in the control hybridization. For 
the copy number analysis, the ratios were normalized on the basis of the 
distribution of ratios of all targets on the array and for the expression analysis 
on the basis of 88 housekeeping genes, which were spotted four times onto the 
array. Low quality measurements {U., copy number data with mean reference 
intensity <I00 Huorcsccnl units, and expression data widi both test and 
reference intensity <I00 fluorescent units and/or with spot size <50 units) 



were excluded from the analysis and were treated as missing values. The 
distributions of fluorescence ratios were used to dcfme outpoints for increased/ 
decreased copy number. Genes with CGH ratio >1.43 (representing the upper 
5% of the CGH ratios across all experiments) were considered to be amplified, 
and genes with ratio <0.73 (representing the lower 5%) were considered to be 
deleted. 

Statistical Analysis of CGH and cDNA Microarray Data, To evaluate 
the influence of copy number alterations on gene expression, we applied the 
following statistical approach. CGH and cDNA calibrated intensity ratios were 
log-transformed and noimalized using median centering of the values m each 
cell line Furthermore. cDNA ratios for each gene across aU 14 cell Imcs were 
median centered. For each gene, the CGH data were represented by a vector 
that was labeled 1 for amplification (ratio. >1 .43) and 0 for no amplification. 
Amplification vws correlated with gene expression using the signaUo-noise 
statistics (1). We calculated a weight, for each gene as follows: 

m,, - m^ 

where m,,. <r,» and m^. denote the means and SDs for the expression 
levels for amplified and nonamplified cdl fmes. respectively. To assess the 
statistical significance of each weight, we performed 10,000 random permu- 
tations of the label vector. The probability diat a gene had a larger or equal 
weight by random permutation than the original weight was denoted by a. A 
low a (<0.05) indicates a strong association between gene expression and 

amplification. . 

Genomic LocaUzation of cDNA Clones and Amplicon Mappmg. Each 
cDNA clone on the micrcianay was alssigned to a Unigene cluster using the 
Unigene BuUd 141 * A database of genomic sequence alignment mformauon 
for mRNA sequences was created from the August 2001 freeze of the Uni- 
versity of aiifomia Santa Crai's GoldenPath database.' The chromosome and 
bp positions for each cDNA clone were then retrieved by relating these data 
sets Amplicons were defined as a CGH copy number ratio >2.0 in at least two 
adjacent clones in two or more cell lines or a CGH ratio >2.0 in at least three 
adjacent clones in a single cell line. The amplicon start and end positions were 



* Internet address: bap-y/rcscarcb.nhgrLnfl».gOY/rmciMiTay/dowiJo»dablc_c^^^ 
' Internet address: www.gciiomc.ucsc.edu. 
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GENE EXPRESSION PATTERNS W BREAST CANCER 

CGH were validated, with lq2U nql2-q21.2, 17q22-q23. 20ql3,l, 
and 20ql3.2 regions being most commonly amplified. Furthemiore, 
the boundaries of these ampHcons were precisely delineated. In ad- 
diUon. novel amplicons were identified at 9pl3 (38.6«9,25 Mb), 
and 17q21.3 (52.47-55.80 Mb). 

Direct Identification of Putative Amplification Target Genes. 
The cDNA/CGH microarray technique enables the direct correla- 
tion of copy number and expression data on a gene-by-gene basis 
throughout the genome. We directly annotated high-rcsolution 
CGH plots with gene expression data using color coding. Fig, 2C 
shows that most of the amplified genes in the MCF-7 breast cancer 
cell line at lpl3. I7q22-q23. and .20ql3 were highly overex- 
pressed, A view of chromosome 7 in the MDA.468 cell hne 
implicates EGFR as the most highly overexpressed and amplified 
gene at 7pl l-pl2 (Fig. 3^). In BT-474, the two known amplicons 
at 17ql2 and 17q22-q23 contained numerous highly overex- 
pressed genes (Fig. 35). In addition, several genes, including the 
homeobox genes H0XB2 and HOXB7, were highly amplified in a 
previously undescribed independent araplicon at I7q21.3. H0XB7 
was systematically amplified (as validated by FISH, Fig. 3BJnset) 
as well as overexpressed (as verified by RT-PCR, data not shown) 
in BT-474 UACC812, and ZR-75-30 cells. Furthermore, this novel 
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extended to include nci^boring nonamplificd clones (ratio. <1.5), The ara- 
plicon size deterauoation was partially dependent on local clone density. 

FISH, Dual-color interphase FISH to breast cancer cell Iine$ was done as 
described (17). Bacterial artificial chromosome clone RP1I-361K8 was la- 
beled with SpcctmmOrange (Vysis. Downers Grove, IL), and Spcctmm- 
Orangclabcled probe for EGFR was obtained from Vysis. SpectrumGrccn- 
labeled chromosome 7 and 17 centromere probes (Vysis) were used as a 
reference. A ussue microarray containing 612 formalin-fixed, parafTin-embcd- 
ded primary breast cancers (17) was applied in FISH analyses as described 
(18) The use of these specimens was approved by the Ethics Commmce of a»c 
University of Basel and by the NIH. Specimens conuining a 2-fold or higher 
increase in the number of test probe signals, as compared with correspondir^ 
centromere signals, in at least 10% of the tumor cells were considered to be 
amplified. Survival analysis was performed using the Kaplan-Mcicr method 
and the log-rank test. . i • 

RT-PCR. The H0XB7 expression level was detennined relative to 
CAPDH, Reverse transcription and PCR amplification were perfomied usmg 
Access RT-PCR System (Promega Corp.. Madison. WI) wiA 10 ^^f^^^ 
as a template. H0XB7 primers were i'-GAGCAGAGGGACTCGGACTT-S 
and 5'-GCGTCAGGTAGCGATTGTAG-3'. 



RESULTS 

Global Effect of Copy Number on Gene Expression. 13,824 
arrayed cDNA clones were applied for analysis of gene expression 
and gene copy number (CGH microanays) in 14 breast cancer cell 
Unes. The results illustrate a considerable influence of copy number 
on gene expression patterns. Up to 44% of the highly amplified 
transcripts (CGH ratio, >2.5) were overexpressed (ie.. belonged to 
the global upper 7% of expression ratios), compared with only 6% for 
genes with normal copy number levels (Fig. 1 A). Conversely, 10.5% 
of the transcripts with high-level expression (cDNA rauo, >10) 
showed increased copy number (Fig. IB). Low-level copy number 
increases and decreases were also associated with similar, although 
less dramatic, outcomes on gene expression (Fig. 1). 

Identification of Distinct Breast Cancer Amplicons. Base-pair 
locations obtained for 1 1,994 cDNAs (86.8%) were used to plot copy 
number changes as a funcrion of genomic position (Fig. 2. Supple- 
ment Fig. A), The average spacing of clones throughout the genome 
was 267 kb. This high-resolution mapping identified 24 independent 
breast cancer amplicons, spanning from 0,2 to 12 Mb of DNA (Table 
I). Several amplification sites detected previously by chromosomal 
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Fig. 4, List of 50 genes with s stotisUcally 
significtttt correlation (o value <a05) between 
gene copy number and gene expression. Name, 
chromosomal location, and the a value for each 
gene ait indicated. The genes have been ordered 
accordtng to their position in the genome. The color 
maps on the right illustrate the copy number and 
expression ratio patterns in the 14 cell lines. The 
key to the color code is riiown at the bottom of the 
graph. Gray squares, missing values. The complete 
list of 270 getics is shown in supplemental Fig. B. 
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amplification was validated to be present in 10,2% of 363 primary 
breast cancers by FISH to a tissue microarray and was associated 
with poor prognosis of the patients (P = 0.001). 

Statistical Identincation and Characterization of 270 Highly 
Expressed Genes In Araplicons. Statistical comparison of expres- 
sion levels of all genes as a function of gene amplification identified 
270 genes whose expression was significantly influenced by copy 
number across all 14 cell lines (Fig. 4, Supplemental Fig. B). Accord- 
ing to the gene ontology data,* 91 of the 270 genes represented 
hypothetical proteins or genes with no functional annotation, whereas 
179 had associated functional information available. Of these, 151 
(84%) arc implicated in apoptosis. cell proliferation, signal transduc- 
tion, and transcription, whereas 28 (16%) had functional annotations 
that could not be directly linked with cancer. 



DISCUSSION 

The importance of recurrent gene and chromosome copy number 
changes in the development and progression of solid tumors has been 
characterized in >1000 publications applying CGH' (9, 10), as well 
as in a large number of other molecular cytogenetic, cytogenetic, and 
molecular genetic studies. The effects of these somatic genetic 
changes on gene expression levels have remained largely unknown, 
although a few studies have explored gene expression changes occur- 
ring in specific ampHcons (15, 19-21). Here, we applied genome- 
wide cDNA microarrays to identify transcripts whose expression 
changes were attributable to underlying gene copy number alterations 
in breast cancer. 

The overall impact of copy number on gene expression patterns was 
substantial with the most dramatic effects seen in the case of high- 



Imcmct ad<lrcss: http://www.gcneoniology.ofg/. 



• Internet address: http://www.ncbi.nlrojjai.gov/cntrcz. 
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level copy number increase. Low-level copy number gams and losses 
also had a significant influence on expression levels of genes m the 
regions affected, but these effects were more subtle on a gene-by-gene 
b^is than those of high-level amplifications. However, the impact of 
low-lcvel gains on the dysregulalion of gene expression pattems m 
cancer may be equally important if not more important than that of 
hioh-level amplifications. Aneuploidy and low-level gams and losses 
of%hromosomal arms represent the most common types of genetic 
alterations m breast and other cancers and, therefore, have an mflu- 
ence on many.genes. Our results in breast cancer extend the recent 
studies on the impact of aneuploidy on global gene expression pat- 
terns in yeast cells, acute myeloid leukemia, and a prostate cancer 
model system (22-24). ^ . t, 

The CGH raicroarray analysis identified 24 mdepcndcnt breast 
cancer amplicons. We defined the precise boundaries for many^am- 
plicons detected previously by chromosomal CGH (9. 10, 25, 26) and 
also discovered novel amplicons that had not been detected previ- 
ously, presumably because of their small size (only 1-2 Mb) or close 
proximity to other larger amplicons. One of these novel amphcons 
involved the homeobox gene region at 17q21.3 and led to the over- 
expression of x\i^ HOXB7 and H0XB2 genes. The homeodomam 
transcription factors are known to be key regulators of embryomc 
aevelopment and have been occasionally reported to undergo aberrant 
expression in cancer (21, 28). H0XB7 transfection induced cell pr<>. 
liferation in melanoma, breast, and ovarian cancer cells and mcrcased 
tumorigenicity and angiogencsis in breast cancer (29-32). The pres- 
ent results imply that gene amplification may be a prominent mech- 
anism for overexpressing H0XB7 in breast cancer and suggest that 
H0XB7 contributes to tumor progression and confers an aggressive 
disease phenotype in breast cancer. This view is supported by our 
finding of amplification of in 10% of 363 primary breast 
cancers, as well as an association of amplification with poor prognosis 
of the patients. 

We carried out a systematic search to identify genes whose 
expression levels across all 14 cell lines were attributable to 
amplification status. Statistical analysis revealed 270 such genes 
(representing -2% of all genes on the array), including not only 
previously described amplified genes, such as HER-2, MYC, 
EGFR, ribosomal protein s6 kinase^ and AIB3, but also numerous 
novel genes such as NRAS-related gene {1^^^). syndecan-l (8q22). 
and bone morphogenic protein (20ql3,l). whose activation by 
amplification may similarly promote breast cancer progression. 
Most of the 270 genes have not been implicated previously m 
breast cancer development and suggest novel pathogenetic inech- 
anisms. Although we would not expect all of them to be causally 
involved, it is intriguing that 84% of the genes with associated 
functional information were implicated in apoptosis. cell prolifer- 
ation signal transduction, transcription, or other cellular processes 
that could directly imply a possible role in cancer progression. 
Therefore, a deuiled characterization of these genes may provide 
biological* insights to breast cancer progression and might lead to 
ihe development of novel therapeutic strategics. 

In summary, we. demonstrate application of cDNA microarrays 
to the analysis of both copy number and expression levels of over 
12 000 transcripts throughout the breast cancer genome, roughly 
once every 267 kb. This analysis provided: (a) evidence of a 
. prominent global influence of copy number changes on gene 
expression levels; (6) a high-resolution map of 24 "^dependent 
amplicons in breast cancer, and (c) identification of a set of 270 
genes, the overexpression of which was statistically attributable to 
gene amplification. Characterization of a novel amplicon at 
17q2h3 implicated amplification and overexpression of the 
H0XB7 cene in breast cancer, including a clinical association 
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between HOXB7 amplification and poor patient prognosis. Overall, 
our results illustrate how the identification of genes activated by 
oene amplification provides a powerful approach to highlight 
genes with an important role in cancer as well as to pnoriuze and 
validate putative targets for therapy development. 
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Microarray analysis reveals a major direct role of 
DNA copy number alteration in the transcriptional 
program of human breast tumors 
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Robert Tibshiraniss. pavid Botsteir»l. Anne-Use Barresen-Dale*. and Patrick O. Brown 

Comprehensive Cancer Center. University of North Carolina. Chapel H.ll, NC 27599 



Contributed by Patrick O. Brown. August 6. 2002 
Genomic DNA copy number alterations are key genetic events in 
the development and progression of human cancers. Here we 
report a genome-wide microarray comparative genomic hybrid- 
ization (array CGH) analysis of DNA copy number v^»r?ation m 
a series of primary human breast tumors. We have profiled DNA 
copy number alteration aaoss 6.691 mapped human genes, in 44 
predominantly advanced, primary breast tumors and 10 breast 
cancer cell lines. While the overall patterns of DNA amplification 
and deletion corroborate previous cytogenetic studies, the high- 
resolution (gene-by-gene) mapping of amplicon boundaries and 
the quantitative analysis of amplicon shape provide significant 
improvement in the localization of candidate oncogenes. ParaHel 
microarray measurements of mRNA levels reveal the remarkable 
degree to which variation in gene copy number contributes to 
variation in gene expression in tumor cells. Specifically, we find 
that 62% of highly amplified genes show moderately or highly 
elevated expression, that DNA copy number Influences gene ex- 
pression across a wide range of DNA copy number alterations 
(deletion, low-, mid- and high-level amplification), that on average, 
a 2-foId change in DNA copy number Is associated with a corre- 
sponding 1.S-f old change in mRNA levels, and that overall, at least 
12% of all the variation in gene expression among the breast 
tumors is directly attributable to underiying variation in gene copy 
number. These findings provide evidence that widespread DNA 
copy number alteration can lead directly to global deregulation of 
gene expression, which may contribute to the development or 
progression of cancer. 

Conventional cytogenetic techniques, including comparative 
genomic hybridization (CGH) (1). have led to the identifi- 
cation of a number of recurrent regions of DNA copy iiumbcr 
alteration in breast cancer cell lines and tumors (2-4). While 
some of these regions contain known or candidate oncogenes 
[e.I FGFRl (8pll), MYC (8q24), CCNDl (nql3), ERBB2 
07ql2). and ZNF2X7 (20ql3)l and tumor suppressor genes 
fRBl {13ql4) and TP53 (17pl3)], the relevant gene(s) withm 
other regions (e.g.. gain of Iq. 8q22. and 17q22-24. and loss of 
8p) remain to be identified. A high-resolution genome-wide 
map, delineating the boundaries of DNA copy number alter- 
ations in tumors, should facilitate the localization and identifi- 
cation of oncogenes and tumor suppressor genes m breast 
cancer. In this study, we have created such a map, usmg 
array-based CGH (5-7) to profile DNA copy number alteraUon 
in a series of breast cancer cell lines and primary tumors. 

An unresolved question is the extent to which the widespread 
DNA copy number changes that we and others have idcntrfied 
in breast tumors alter expression of genes within involved 
regions. Because we had measured mRNA levels in parallel iq 
the same samples (8), using the same DNA microarrays, we had 
an opportunity to explore on a genomic scale the relationship 
between DNA copy number changes and gene expression. From 

www4>nas.org/<gi/doi/l0.l073/pnas-16247l999 



this analysis, we have identified a significant impact of wide- 
spread DNA copy number alteration on the transcriptional 
programs of breast tumors. 

Materia and Methods - 
Tumors and Cell Line*. Primary breast tumors were predomm^tly 
large (>3 cm), intermediate-grade, infiltrating ductal carano- 
m^S. with more than 50% being lymph "Ode positive T^e 
fraction of tumor cells within specimens averaged at least 50%, 
Details of individual tumors have been pubhshed (H. V). ana 
are summarized in Table 1, which is published ^ supporting 
information on the PNAS web «te. jvww.pnas.org. Bi^^^^^ 
cell lines were obtained from the American Type Culture 
Collection. Genomic DNA was isolated cither usmg Qiagen 
genomic DNA columns, or by phenol/chloroform extraction 
followed by ethanol precipitation. 

DNA Ubeling and Microarray Hybridizations. Genomic DNA label- 
ing and hybridizations were performed essentially as described 
in Pollack er aL (7). with slight modifications. Two micrograms 
of DNA was labeled in a total volume of 50 microliters a^^^^ 
volumes of all reagents were adjusted accordingly. L)NA 
(from tumors and cell lines) was f luorescently labeled (Cy5) and 
hybridized to a human cDNA microarray containing^ 6,691 
different mapped human genes (i.e.. UniGene clusters). The 
"reference" (labeled with Cy3) for each hybridi^tion was nor- 
mal female leukocyte DNA from a single donor The fabrication 
of cDNA microarrays and the labeling and hybridization of 
mRNA samples have been described (8). 

Data Analysis and Map Positions, Hybridized arrays were scanned 
on a GenePix scanner (Axon Instruments. Foster City, CA). and 
fluorescence ratios (test/reference) calculated using scan aly2:e 
software (available at http://rana.lbl.gov). Fluorescence ratios 
were norinalized for each array by setting the average log 
fluorescence ratio for all array elements equa^ to 0- ^^5^"^^- 
mcnts with fluorescence intensities more than 20% above back- 
ground were considered reliable, DNA copy number profiles 
ihzi deviated significandy from background ratios measured m 
normal genomic DNA control hybridizations were mtcfPf 
evidence of real DNA copy number alteration (see Estwianng 
Significance of Altered Fluorescence Ratios m the supporting 
information). When indicated. DNA copy number profiles are 
displayed as a moving average (symmetnc 5-nearest "^'ghbors^ 
Map positions for arrayed human cDNAs were assigned by 



Abbreviation: CGH. compantlve genomic hybridHatton. 
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identifying the starting position of the best and longest match of 
any DN A sequence represented in the corresponding UniGene 
cluster (10) against the **GoIden Path" genome assembly 
(http://gcnomc.ucsc.edu/; Oct 7. 2000 Freeze). For UniGcne 
clusters represented by multiple arrayed elements, mean fluo- 
rescence ratios (for all elements representing the same UniGene 
cluster) are reported. For mRNA measurements, fluorescence 
ratios are "mean-centered" (i.e., reported relative to the mean 
ratio across the 44 tumor samples). The data set descnbed here 
can be accessed in its entirety in the supporting information. 

Results 

We performed CGH on 44 predominantly locally advanced, 
primary breast tumors and 10 breast cancer cell lines, usmg 
cDNA microarrays containing 6.691 different mapped human 
genes (Fig. la; also see Materials and Methods for details of 
microarray hybridizations). To take full advantage of the im- 
proved spatial resolution of array CGH, we ordered (fluores- 
cence ratios for) the 6,691 cDNAs according to the "Golden 
Path" (http://genomc.ucsc.edu/) genome assembly of the draft 
human genome sequences (11). In so doing, arrayed cDNAs not 
only themscWes represent genes of potential interest (e.g., 
candidate oncogenes within amplicons), but also provide precise 
genetic landmarks for chromosomal regions of amplification and 

I ¥tfww.po«.Ofg/cgI/doi/10.1073/pna$. 162471999 



deletion. Parallel analysis of DNA from cell lines contammg 
different numbers of X chromosomes (Fig. 16), as we did before 
m demonstrated the sensitivity of our method to detect smgle- 
copy loss (45, XO), and 1^- (47,XXX). 2- (48,XXXX) or 
2.5-fold (49,XXXXX) gains (also sec Fig. 5, which is published 
as supporting information on the PNAS web site). Fluorescence 
ratios were linearly proportional to copy number ratios, which 
were slightly underestimated. In agreement with previous ob- 
servations (7). Numerous DNA copy number alterations were 
evident in both the breast cancer cell lines and pnmary turaore 
(Fig. la), detected in the tumors despite the presence of euploid 
non-tumor cell types; the magnitudes of the observed changes 
were generally lower in the tumor samples. DNA copy-number 
alterations were found in every cancer cell line and tumor, and 
on every human chromosome in at least one sample. Recurrent 
regions of DNA copy number gain and loss were readily iden- 
tifiable. For example, gains within Iq, 8q, 17q, and 20q were 
observed in a high proportion of breast cancer cell lines/tumors 
(90%/69%, 100%/47%, 100%/60%, and 90%/44% respecttw- 
W), as were losses within Ip. 3p, 8p, and 13q (80%/24%, 
80%/22%, 80%/22%, and 10%/IS%, respectively), consistent 
with published cytogenetic studies (refe. 2-4; a complete hsUng 
of gains/losses is provided in Tables 2 and 3, which are published 
as supporting information on the PNAS vfcb site). The total 
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number of genomic alterations (gains and losses) was found to 
be significantly higher in breast tumors that were high grade {P - 
0.008). consistent with published GGH data (3), estrogen recep- 
tor negative (P - 0.04). and harboring TP53 mutations {P = 
0.0006) (sec Table 4, which is published as supporting informa- 
tion on the PNAS web site). . 

The improved spatial resolution of our array CGH analysis is 
illuslraled for chromosome 8, which displayed extensive DNA 
copy number alteration in our series. A detailed view of the 
variation in the copy number of 241 genes mapping to diromo- 
some 8 revealed multiple regions of recurrent amphficaiion; 
each of these potentially hartwrs a different known or previously 
uncharacterizcd oncogene (Fig. 2a). The complexity of amphcon 
structure is most easily appreciated in the breast cancer cell Imc 
SKBR3. Although a conventional CGH analysis of 8q in SKBR3 
identified only two distinct regions of amplification (12), we 
observed three distinct regions of high-level amplification (la- 
beled 1-3 in Fig. 2b). For each of these regions we can define the 
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boundaries of the interval recurrently amplified in the tumors we 
examined; in each case, known or plausible candidate oncogenes 
can be identified (a description of these regions, as well as the 
recunently amplified regions on chromosomes 17 and 20. can be 
found in Figs. 6 and 7. which are published as supportmg 
information on the PNAS web site). /.Am 

For a subset of breast cancer cell lines and tumcirs and 37, 
respectively), and a subset of arrayed genes (6,095), 
levels were quantitatively measured in parallel by usmg cDNA 
microarrays (8). The parallel assessment of mRNA levels is 
usehil in the interpretation of DNA copy number changes. For 
example, the highly amplified genes that are also highly ex- 
pressed arc the strongest candidate oncogenes within an ampli- 
con. Perhaps more significantly, our parallel analysis of DNA 
copy number changes and mRNA levels provides ;s the oppor- 
tunity to assess the global impact of widespread DNA copy 
number alteration on gene expression in tumor cells. 

A strong influence of DNA copy number on gene expression 
is evident in an examination of the pseudocolor representations 
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of DNA copy number and mRNA levels for genes on chromo- 
some 17 (Fig. 3). The overall patterns of gene amplification and 
elevated gene expression arc quite concordant; i.e., a significant 
fraction of highly amplified genes appear to be correspondingly 
highly expressed. The concordance between high-level ampli- 
cation and increased gene expression is not restricted to chro- 
mosome 17. Genomc-wide, of 117 high-level DNA araplifica. 
tions (fluorescence ratios >4. and reprcsentmg 91 different 
genes), 62% (representing 54 different genes; see Table 5, which 
is published as supporting information on the PNAS web sue) 
are found assodated with at least moderately elevated mRNA 
levels (mean-centered fluorescence ratios >2), and 42% (rep- 
resenting 36 different genes) are found associated with compa- 
rably highly elevated mRNA levels (mcan^ntered fluorescence 

'^To^ determine the extent to which DNA deletion and lower- 
level amplification (in addition to high-level amplification) are 
also associated with corresponding alterations in mRNA levels, 
we performed three separate ana^ on the complete data set 
(4 cell lines and 37 tumors, across 6.095 genes). First, we 
determined the average mRNA levels for eadi of five classes 
of genes, representing DNA deletion, no change, and low-, 
medium-, and high-level amplification (Fig. 4a). For both the 

12966 I www.pfMi$.org/cgi/doi/10,i073/pna$.1 62471999 



breast cancer cell lines and tumors, average mRNA levels 
tracked with DNA copy number across all five classes, m a 
statistically significant fashion (P values ^^J^fj^^if Student s 
rtestscomparingadjacentclasscs:(»^^n^^^^^ Vv 10-^^' 

5 X 10-^ 1 X 10-^; tumors, 1 x 10 « 1 x 10 " . 5 x 10 , 
1 X 10-*). A linear regression of the average logP^A copy 
number), for each class, against avcra^ ^""^^^ma J^v 
demonstrated that on average, a 2-fold change m DNA copy 
number was accompanied by 1.4- and 1.5-fold changes m mRNA • 
level for the breast cancer cell lines and tumors, respectively (Fig. 
4a regression line not shown). Second, wc characterized the 
distribution of the 6.095 correlations between DNA copy num- 
ber and mRNA level, each across the 37 tumor samples (Fig. 4d). 
The distribution of correlations forms a normal-shaped curve, 
but with the peak markedly shifted in the posiUye direction from 
zero. This shift is statistically significant, as evidenced in a plot 
of observed vs. expected correlations (Fig. 4c). and reflects a 
pervasive global influence of DNA copy nuniber alteraUons on 
gene expr^sion. Notably, the highest correlations between DN^^ 
copy number and mRNA level (the right taU of the d>stnbut|on 
in Fir 46) comprise both amplified and deleted genes (data not 
shown). Third, we used a linear regression model to estimate the 
fraction of all variation measured in mRNA levels among the 37 
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tumors that could be attributed to underlying variation in DNA 
copy number. From this analysis, we estimate that, overall, about 
7% of all of the observed variation in mRNA levels can be 
explained directly by variation in copy number of the altered 
genes (Fig. 4d), We can reduce the effects of experimental 
measurement error on this estimate by using only that fraction 
of the data most reliably measured (fluorescence intensity/ 
background >3); using that data, our estimate of the percent 
variation in mRNA levels directly attributed to variation in gene 
copy number increases to 12% (Fig. 4^0- This still undoubtedly 
represents a significant underestimate, as the observed variation 
in global gene expression is affected not only by true variation in 
the expression programs of the tumor cells thcmseNes, but also 
by the variable presence of non-tumor cell types within clinical 
samples. 

Discussion 

This genome-wide, array CGH analysis of DNA copy number 
alteration in a scries of human breast tumors demonstrates the 
usefuhiess of defining amplicon boundaries at high resolution 
(gene-by-gene), and quantitatively measuring amplicon shape, to 
assist in locating and identifying candidate oncogenes. By ana- 
lyzing mRNA levels in parallel, we have also discovered that 
changes in DNA copy number have a large, pervasive, direct 
effect on global gene expression patterns in both breast cancer 



cell lines and tumors. Although the DNA microarrays used in our 
analysis may display a bias toward characterized and/or highly 
expressed genes, because we are examining such a large fraction 
of the genome (approximately 20% of all human genes), and 
because, as delated above, we arc likely underestimating the 
contribution of DNA copy number changes to altered gene 
expression, we believe our findings arc likely to be generalizable 
(but would nevertheless still be remarkable if only applicable to 
this set of -6,100 genes). 

In budding yeast, ancuploidy has been shown to result m 
chromosome-wide gene expression biases (13). Two recent 
studies have begun to examine the global relationship between 
DNA copy number and gene expression in cancer cells. In 
agreement with our findings, Phillips ei aL (14) have shown that 
with the acquisition of tumorigcnicity in an immortalized pros- 
tate epithelial cell line, new chrompsomal gains and losses 
resulted in a statbtically significant respective increase and 
decrease in the average expression level of involved genes. In 
contrast, Platzer et aL (15) recently reported that in metastatic 
colon tumors only -4% of genes within amplified regions were 
found more highly (>2-fbld) expressed, when compared with 
normal colonic epithelium. This report differs substantially from 
our finding that 62% of highly amplified genes in breast cancer 
exhibit at least 2-fold increased expression. These contrasting 
findings may reflect methodological differences between the 
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studies. For example, the study of Platzer et al (15) may have 
systematically under-measured gene expression changes. In this 
regard it is remarkable that only 14 transcripts of many thousand 
residing within unamplified chromosomal regions were found to 
exhibit at least 4-fold altered expression in metastatic colon 
cancer. Additionally, their reliance on lower-resolution chromo- 
somal CGH may have, resulted in poorly delimiting the bound- 
aries of high-complexily amplicons, effectively overcalling re- 
gions with amplification. Altcniatively, the contrasting findings 
for amplified genes may represent real biological differences 
between breast and metastatic colon tumors; resolution of this 
issue will require further studies. 

Our finding that widespread DNA copy number alteration has 
a large, pervasive and direct effect on global gene expression 
patterns in breast cancer has several important implications. 
First, this finding supports a high degree of copy number- 
dependent gene expression in tumors. Second, it suggests that 
most genes are not subject to specific autoregulation or dosage 
compensation. Third, this finding cautions that elevated expres- 
sion of an amplified gene carmot alone be considered strong 
independent evidence of a candidate oncogene's role in turaor- 
igenesis. In our study, fully 62% of highly amplified genes 
demonstrated moderately or highly elevated expression. This 
highlights the importance of high-resolution mapping of ampli- 
con boundaries and shape (to identify the "driving" gene(s) 
within amplicons (16)], on a large number of samples, in addition 
to functional studies. Fourth, this finding suggesU that analyzing 
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the genomic distribution of expressed genes, even within existing 
microarray gene expression data sets, may permit the inference 
of DNA copy number aberration, particularly aneuploidy (where 
gene expression can be averaged across large chromosomal 
regions; sec Fig. 3 and supporting information). Fifth, this 
finding implies that a substantial portion of the phenotypic 
uniqueness (and by extension, the heterogeneity in clinical 
behavior) among patients* tumon may be traceable to underly- 
ing variation in DNA copy number. Sixth, this finding supports 
a possible role for widespread DNA copy number alteration in 
tumorigcncsis (17, 18), beyond the amplification of specific 
oncogenes and deletion of spedfic tumor suppressor genes. 
Widespread DNA copy number alteration, and the concomitant 
widespread imbalance in gene expression, might disrupt aitical 
stochiometric relationships in cell metabolism and physiology 
(e.g., proteosome. mitotic spindle), possibly promoting further 
chromosomal instability and directly contributing to tumor 
development or progression. Finally, our findings suggest the 
possibility of cancer therapies that exploit specific or global 
imbalances in gene expression 'in cancer. 
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Each year, over 1 82.000 women in the United States are 
diagnosed with breast cancer, and approximately 45,000 die 
of the diseased Incidence appears to be increasing in the 
United States at a rate of roughly 2% per year. The reasons 
for the increase are unclear, but non-genetic risk factors appear 
to play a large role.2 

Five-year survival rates range from approximately 65%- 
85%, depending on demographic group, with a significant 
percentage of women experiencing recurrence of their cancer 
within 10 years of diagnosis. One of the factors most predic- 
tive for recurrence once a diagnosis of breast cancer has been 
made is the number of axillary lymph nodes to which tumor 
has metastasized. Most node-positive women are given adju- 
vant therapy, which increases their survival. However, 20%- 
30% of patients without axillary node involvement also 
develop recurrent disease, and the difficulty lies in how to iden- 
tify this high-risk subset of patients. These patients could 
benefit from increased surveillance, early intervention, and 
treatment. 

Prognostic markers currently used in breast cancer recur- 
rence prediction include tumor size, histological grade, steroid 
hormone receptor status, DNA ploidy, proliferative index, and 
cathepsin D status. Expression of growth factor receptors and 
over-expression of the HER-2/neu oncogene have also been 
identified as having value regarding treatment regimen and 
prognosis. 

HER-2/neu (also known as c-erbB2) is an oncogene that 
encodes a transmembrane glycoprotein that is homologous 
to, but distinct from, the epidermal growth factor receptor. 
Numerous studies have indicated that high levels of expres- 
sion of this protein are associated with rapid tumor growth, 
certain forms of therapy resistance, and shorter disease-free 
survival. The gene has been shown to be amplified and/or 
overexpressed in 1 0%-30% of invasive breast cancers and in 
40%-60% of intraductal breast carcinoma.^ 

There are two distinct FDA-approved methods by which 
HER-2/neu status can be evaluated: immunohistochemistry 
(IHC, HcrccpTest™) and FISH (fluorescent in situ hybridiza- 
tion, PathVysion""** Kit). Both methods can be performed on 
archived and current specimens. The first method allows visual 
assessment of the amount of HER-2/neu protein present on 
the cell membrane. The latter method allows direct quantifi- 
cation of the level of gene amplification present in the tumor, 
enabling differentiation between low- versus high-amplifica- 
tidh. At least one study has demonstrated a difference in 



recurrence risk in women younger than 40 years of age for 
low- versus high-amplified tumors (54.5% compared to 
85.7%); this is compared to a recurrence rate of 16.7% for 
patients with no HER-2/neu gene amplification.^ HER-2/neu 
status may be particularly important to establish in women with 
small (^1 cm) tumor size. 

The choice of methodology for determination of HER-2/ 
neu status depends in part on the clinical setting. FDA approval 
for the Vysis FISH test was granted based on clinical trials 
involving 1549 node-positive patients. Patients received one 
of three different treatments consisting of different doses of 
cyclophosphamide, Adriamycin, and 5-fluorouracil (CAF). 
The study showed that patients with amplified HER-2/neu 
benefited from treatment with higher doses of adriamycin- 
based therapy, while those with normal HER-2/neu levels did 
not. The study therefore identified a sub-set of women, who 
because they did not benefit from more aggressive treatment, 
did not need to be exposed to the associated side effects. In 
addition, other evidence indicates that H£R-2/neu amplifica- 
tion in node-negative patients can be used as an independent 
prognostic indicator for early recurrence, recurrent disease at 
any time and disease-related death.^ Demonstration of HER- 
2/neu gene amplification by FISH has also been shown to be 
of value in predicting response to cheniotherapy in stage-2 
breast cancer patients. 

Selection of patients for Herceptin® (Trastuzumab) mono- 
clonal antibody therapy, however, is based upon demonstra- 
tion of HER-2yneu protein overexpression using HercepTest™. 
Studies using Herceptin^ in patients with metastatic breast 
cancer show an increase in time to disease progression, 
increased response rate to chemotherapeutic agents and a small 
increase in overall survival rate. The FISH assays have not yet 
been approved for this purpose, and studies looking at response 
to Herceptin® in patients with or without gene amplification 
status determined by FISH are in progress. 

In general, FISH and IHC results correlate well However, 
subsets of tumors are found which show discordant results; 
i.e., protein overexpression without gene amplification or lack 
of protein overexpression with gene amplification. The clini- 
cal significance of such results is unclear. Based on the above 
considerations, HER-2/neu testing at SHMC/PAML will uti- 
lize immunohistochemistry (HercepTest^) as a screen, fol- 
lowed by FISH in IHC-negative cases. Alternatively, either 
method may be ordered individually depending on the clini- 
cal setting or clinician preference. 
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CPT code information 

HER-2/tteu via IHC 

88342 (including interpretive report) 

HER-2/neu via FISH 

88271 x2 Molecular cytogenetics, DNA probe, each 
88274 Molecular cytogenetics, interphase in situ hybrid- 
ization, analyze 25-99 cells 
8829 1 Cytogenetics and molecular cytogenetics, interpre- 
tation and report 



Procedural Information 

Immunohistochemistry is performed using the FDA-approved 
DAKO antibody kit, Herceptest^*. The DAKO kit contains 
reagents required to complete a two-step immunohisto- 
chemical staining procedure for.routinely processed, paraffin- 
embedded specimens. Following incubation with the primary 
rabbit antibody to human HER-2/neu protein, the kit employs 
a ready-to-use dexlran-based visualization reagent. This re- 
agent consists of both secondary goat anti-rabbit antibody 
molecules with horseradish peroxidase molecules linked to a 
common dextran polymer backbone, thus eliminating the need 
for sequential application of link antibody and peroxidase 
conjugated antibody. Enzymatic conversion of the subse- 
quently added chromogen results in formation of visible 
reaction product at the antigen site. The specimen is then coun- 
terstained; a pathologist using light-microscopy interprets 
results. 

FISH analysis at SHMC/PAML is performed usmg the 
FDA-approved PaihVysion™ HER-2/neu DNA probe kit, pro- 
duced by Vysis, Inc. Formalin fixed, parafFm-embcdded breast 
tissue is processed using routine histological methods, and then 
slides are treated to allow hybridization of DNA probes to the 
nuclei present in the tissue section. The Pathvysion™ kit con- 
tains two direct-labeled DNA probes, one specific for the 
alphoid repetitive DNA (CEP 1 7, spectrum orange) present at 
the chromosome 17 centromere and the second for the HER- 
2/neu oncogene located at 1 7q 1 1 .2- 1 2 (spectrum green). Enu- 
meration of the probes allows a ratio of the number of copies 
of chromosome 17 to the number of copies of HER-2/neu to 
be obtained; this enables quantification of low versus high 
amplification levels, and allows an estimate of the percentage 
of cells with HER-2/neu gene amplification. The clinically 
relevant distinction is whether the gene amplification is due 
10 increased gene copy number on the two chromosome 17 
homologues normally present or an increase in the number of 
chromosome 17s in the cells. In the majority of cases, ratio 
equivalents less than 2.0 are indicative of a normal/negative 
result, ratios of 2.1 and over indicate that amplification is 
present and to what degree. Interpretation of this data will be 
performed and reported from the Vysis-certificd Cytogenet- 
ics laboratory at SHMC. 
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ABSTRACT Wnt family members are criticial to many 
developmental processes, and components of the Wnt signal- 
ing pathway have been linked to tumorigenesis in familial and 
sporadic colon carcinomas. Here we report the identification 
of two genes, WISP-l and WISP'2, that are up-regulated in the 
mouse mammary epithelial cell line C57MG transformed by 
Wnt-1, but not by Wnt-4. Together with a third related gene, 
WISP'3, these proteins define a subfamily of the connective 
tissue growth factor family. Two distinct systems demon- 
strated WISP induction to be associated with the expression of 
Wnt-1. These included (i) C571VIG cells infected with a Wnt-1 
retroviral vector or expressing Wnt-1 under the control of a 
tetracyline repressible promoter, and («) Wnt-1 transgenic 
mice. The WISP'l gene was localized to human chromosome 
8q24.1-8q24.3. WISP- J genomic DNA was amplified in colon 
cancer cell lines and in human colon tumors and its RNA 
overexpressed (2- to >30-fold) in 84% of the tumors examined 
compared with patient-matched normal mucosa. WISP-S 
mapped to chromosome 6q22-6q23 and also was overex- 
pressed (4- to > 40-fold) in 63% of the colon tumors analyzed. 
In contrast, WISP'2 mapped to human chromosome 20ql2- 
20ql3 and its DNA was amplified, but RNA expression was 
reduced (2- to >30-fold) in 79% of the tumors. These results 
suggest that the WISP genes may be downstream of Wnt-1 
signaling and that aberrant levels of WISP expression in colon 
cancer may play a role in colon tumorigenesis. 



Wnt-1 is a member of an expanding family of cysteine-rich, 
glycosylated signaling proteins that mediate diverse develop- 
mental processes such as the control of cell proliferation, 
adhesion, cell polarity, and the establishment of cell fates (1, 
2), Wnt-1 originally was identified as an oncogene activated by 
the insertion of mouse mammary tumor virus in virus-induced 
mammary adenocarcinomas (3, 4). Although Wnt-1 is not 
expressed in the normal mammary gland, expression of Wnt-1 
in transgenic mice causes mammary tumors (5). 

In mammalian cells, Wnt family members initiate signaling 
by binding to the seven-transmembrane spanning Frizzled 
receptors and recruiting the cytoplasmic protein Dishevelled 
(Dsh) to the cell membrane (1, 2, 6). Dsh then inhibits the 
kinase activity of the normally constitutively active glycogen 
synthase kinase-3)3 (GSK-3^) resulting in an increase in 
)3-catenin levels. Stabilized )3-catenin interacts with the tran- 
scription factor TCF/Lefl, forming a complex that appears in 
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the nucleus and binds TCF/Lefl target DNA elements to 
activate transcription (7, 8). Other experiments suggest that 
the adenomatous polyposis coli (APC) tumor suppressor gene 
also plays an important role in Wnt signaling by regulating 
j3-catenin levels (9). APC is phosphorylated by GSK-3ft binds 
to )3-catenin, and facilitates its degradation. Mutations in 
either APC or ^-catenin have been associated with colon 
carcinomas and melanomas, suggesting these mutations con- 
tribute to the development of these types of cancer, implicating 
the Wnt pathway in tumorigenesis (1). 

Although much has been learned about the Wnt signaling 
pathway over the past several years, only a few of the tran- 
scriptionally activated downstream components activated by 
Wnt have been characterized. Those that have been described 
cannot account for all of the diverse functions attributed to 
Wnt signaling. Among the candidate Wnt target genes are 
those encoding the nodal-related 3 gene, ATirJ, a member of 
the transforming growth factor (TGF)-^ superfamily, and the 
homeobox genes, engrailed, goosecoid, twin (Xtwn), and siamois 
(2). A recent report also identifies c-myc zs a target gene of the 
Wnt signaling pathway (10). 

To identify additional downstream genes in the Wnt signal- 
ing pathway that are relevant to the transformed cell pheno- 
type, we used a PCR-based cDNA subtraction strategy, sup- 
pression subtractive hybridization (SSH) (11), using RNA 
isolated from C57MG mouse mammary epithelial cells and 
C57MG cells stably transformed by a Wnt-1 retrovirus. Over- 
expression of Wnt-1 in this cell line is sufficient to induce a 
partially transformed phenotype, characterized by elongated 
and refractile cells that lose contact inhibition and form a 
multilayered array (12, 13). We reasoned that genes differen- 
tially expressed between these two cell lines might contribute 
to the transformed phenotype. 

In this paper, we describe the cloning and characterization 
of two genes up-regulated in Wnt-1 transformed cells, WISP-l 
and WlSP-2, and a third related gene, WISPS, The WISP genes 
are members of the CCN family of growth factors, which 
includes connective tissue growth factor (CTGF), Cyr61, and 
nov, a family not previously linked to Wnt signaling. 

MATERIALS AND METHODS 

SSH. SSH was performed by using the PCR-Select cDNA 
Subtraction Kit (CLONTECH). Tester double-stranded 

Abbreviations: TGF, transforming growth factor; CTGF, connective 
tissue growth factor; SSH, suppression subtractive hybridization; 
VWC, von Willebrand factor type C module. 

Data deposition: The sequences reported in this paper have been 
deposited in the Genbank database (accession nos. AF100777, 
AF100778, AF100779, AF100780, and AF100781). 
tTo whom reprint requests should be addressed. e*mail: diane@gene. 
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cDNA was synthesized from 2 /xg of polyCA)"^ RNA isolated 
from the C57MG/Wnt-1 cell line and driver cDNA from 2 /ig 
of poly(A)-^ RNA from the parent C57MG cells. The sub- 
tracted cDNA library was subcloned into a pGEM-T vector for 
further analysis. 

cDNA Library Screening. Clones encoding full-length 
mouse WlSP'l were isolated by screening a AgtlO mouse 
embryo cDNA library (CLONTECH) with a 70-bp probe from 
the original partial clone 568 sequence corresponding to amino 
acids 128-169. Clones encoding full-length human WISP-l 
were isolated by screening AgtlO lung and fetal kidney cDNA 
libraries with the same probe at low stringency. Clones en- 
coding full-length mouse and human WlSP-2 were isolated by 
screening a C57MG/Wnt-1 or human fetal lung cDNA library 
with a probe corresponding to nucleotides 1463-1512. Full- 
length cDNAs encoding WISP-3 were cloned from human 
bone marrow and fetal kidney libraries. 

Expression of Human WISP RNA, PCR amplification of 
first-strand cDNA was performed with human Multiple Tissue 
cDNA panels (CLONTECH) and 300 /iM of each dNTP at 
94°C for 1 sec, 62°C for 30 sec, 7T'C for 1 min, for 22-32 cycles. 
WISP and glyceraldehyde-3-phosphate dehydrogenase primer 
sequences are available on request. 

In Situ Hybridization. ^^P-labeled sense and antisense ribo- 
probes were transcribed from an 897-bp PCR product corre- 
sponding to nucleotides 601-1440 of mouse WISP-l or a 
294-bp PCR product corresponding to nucleotides 82-375 of 
mouse WISP'2. All tissues were processed as described (40). 

Radiation Hybrid Mapping. Genomic DNA from each 
hybrid in the Stanford G3 and Genebridge4 Radiation Hybrid 
Panels (Research Genetics, Huntsville, A L) and human and 
hamster control DNAs were PCR-amplified, and the results 
were submitted to the Stanford or Massachusetts Institute of 
Technology web servers. 

Cell Lines, Tumors, and Mucosa Specimens. Tissue speci- 
mens were obtained from the Department of Pathology (Uni- 
versity of Pittsburgh) for patients undergoing colon resection 
and from the University of Leeds, United Kingdom. Genomic 
DNA was isolated (Qiagen) from the pooled blood of 10 
normal human donors, surgical specimens, and the following 
ATCC human cell lines: SW480, CGLO 320DM, HT-29, 
WiDr, and SW403 (colon adenocarcinomas), SW620 (lymph 
node metastasis, colon adenocarcinoma), HCT 116 (colon 
carcinoma), SK-CO-1 (colon adenocarcinoma, ascites), and 
HM7 (a variant of ATCC colon adenocarcinoma cell line LS 
174T). DNA concentration was determined by using Hoechst 
dye 33258 intercalation f luorimetry. Total RNA was prepared 
by homogenization in 7 M GuSCN followed by centrifugation 
over CsCl cushions or prepared by using RNAzol. 

Gene AmpliHcation and RNA Expression Analysis. Relative 
gene amplification and RNA expression of WISPs and c-myc in 
the cell lines, colorectal tumors, and normal mucosa were 
determined by quantitative PCR. Gene-specific primers and 
fluorogenic probes (sequences available on request) were 
designed and used to amplify and quantitate the genes. The 
relative gene copy number was derived by using the formula 
2(^«) where ACt represents the difference in amplification 
cycles required to detect the WISP genes in peripheral blood 
lymphocyte DNA compared with colon tumor DNA or colon 
tumor RNA compared with normal mucosal RNA. The 
d-method was used for calculation of the SE of the gene copy 
number or RNA expression level. The W^/5F-specific signal was 
normalized to that of the glyceraldehyde-3-phosphate dehy- 
drogenase housekeeping gene. All TaqMan assay reagents 
were obtained from Perkin-Elmer Applied Biosystems. 

RESULTS 

Isolation of WISP-l and WISP-2 by SSH. To identify Wnt- 
1-inducible genes, we used the technique of SSH using the 



mouse mammary epithelial cell line C57MG and C57MG cells 
that stably express Wnt-1 (11). Candidate differentially ex- 
pressed cDNAs (1,384 total) were sequenced. Thirty-nine 
percent of the sequences matched known genes or homo- 
logues, 32% matched expressed sequence tags, and 29% had 
no match. To confirm that the transcript was differentially 
expressed, semiquantitative reverse transcription-PCR and 
Northern analysis were performed by using mRNA from the 
C57MG and C57MG/Wnt-1 cells. 

Two of the cDNAs, WISP-l and WISP-2, were differentially 
expressed, being induced in the C57MG/Wnt-1 cell line, but 
not in the parent C57MG cells or C57MG cells overexpressing 
Wnt-4 (Fig. 1 A and 5). Wnt-4, unlike Wnt-1, does not induce 
the morphological transformation of C57MG cells and has no 
effect on jS-catenin levels (13, 14). Expression of WISP-l was 
up-regulated approximately 3-fold in the C57MG/Wnt-1 cell 
line and WISP-2 by approximately 5-fold by both Northern 
analysis and reverse transcription-PCR. 

An independent, but similar, system was used to examine 
WISP expression after Wnt-1 induction. C57MG cells express- 
ing the Wnt-1 gene under the control of a tetracycline- 
repressible promoter produce low amounts of Wnt-1 in the 
repressed state but show a strong induction of Wnt-1 mRNA 
and protein within 24 hr after tetracycline removal (8). The 
levels of Wnt-1 and WISP RNA isolated from these cells at 
various times after tetracycline removal were assessed by 
quantitative PCR. Strong induction of Wnt-1 mRNA was seen 
as early as 10 hr after tetracycline removal. Induction of WISP 
mRNA (2- to 6-fold) was seen at 48 and 72 hr (data not shown). 
These data support our previous observations that show that 
WISP induction is correlated with Wnt-1 expression. Because 
the induction is slow, occurring after approximately 48 hr, the 
induction of WISPs may be an indirect response to Wnt-1 
signaling. 

cDNA clones of human WISP-l were isolated and the 
sequence compared with mouse WISP-L The cDNA sequences 
of mouse and human WISP-l were 1,766 and 2,830 bp in length, 
respectively, and encode proteins of 367 aa, with predicted 
relative molecular masses of •^40,000 (Afr 40 K). Both have 
hydrophobic N-terminal signal sequences, 38 conserved cys- 
teine residues, and four potential N-linked glycosylation sites 
and are 84% identical (Fig. Z4). 

Full-length cDNA clones of mouse and human WISP-2 were 
1,734 and 1,293 bp in length, respectively, and encode proteins 
of 251 and 250 aa, respectively, with predicted relative molec- 
ular masses of 27,000 (Afr 27 K) (Fig. 25). Mouse and human 
WISP-2 are 73% identical. Human WISP-2 has no potential 
N-linked glycosylation sites, and mouse WISP-2 has one at 
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Fig. 1 . WlSP-1 and WISP'2 are induced by Wnt-1, but not Wnt-4, 
expression in C57MG cells. Northern analysis of WISP-l {A) and 
WISP-2 (B) expression in C57MG, C57MG/Wnt-1, and C57MG/ 
Wnt-4 cells. Poly(A)+ RNA (2 tx$) was subjected to Northern blot 
analysis and hybridized with a 70-bp mouse W5P-7-specific probe 
(amino acids 278-300) or a 190-bp W75P-2-specific probe (nucleotides 
1438-1 627) in the 3' untranslated region. Blots were rehybridized with 
human /3-actin probe. 
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Fig. 2. Encoded amino acid sequence alignment of mouse and 
human WISP-1 {A) and mouse and human WISP'2 (B). The potential 
signal sequence, insulin-like growth factor-binding protein (IGF-BP), 
VWC, thrombospondin (TSP), and C-terminal (CT) domains are 
underlined. 

position 197. WISP'2 has 28 cysteine residues that are con- 
served among the 38 cysteines found in WISP-L 

Identification of WJSP-S. To search for related proteins, we 
screened expressed sequence tag (EST) databases with the 
WISP-1 protein sequence and identified several ESTs as 
potentially related sequences. We identified a homologous 
protein that we have called WISP-3. A full-length human 
WISP-3 cDNA of 1,371 bp was isolated corresponding to those 
ESTs that encode a 354-aa protein with a predicted molecular 
mass of 39,293. WISP-3 has two potential N-linked glycosyl- 
ation sites and 36 cysteine residues. An alignment of the three 
human WISP proteins shows that WISP-1 and WISP-3 are the 
most similar (42% identity), whereas WlSP-2 has 37% identity 
with WISP-1 and 32% identity with WISP-3 (Fig. 14). 

mSPs Are Homologous to the CTGF Family of Proteins. 
Human WISP-l, WISP'2, and WISP-3 are novel sequences; 
however, mouse WISP-1 is the same as the recently identified 
Elml gene. Elml is expressed in low, but not high, metastatic 
mouse melanoma cells, and suppresses the in vivo growth and 
metastatic potential of K-1735 mouse melanoma cells (15). 
Human and mouse WISP'2 are homologous to the recently 
described rat gene, rCop-1 (16). Significant homology (36- 
44%) was seen to the CCN family of growth factors. This family 
includes three members, CTGF, Cyr61, and the protoonco- 
gene nov. CTGF is a chemotactic and mitogenic factor for 
fibroblasts that is implicated in wound healing and fibrotic 
disorders and is induced by TGF-^ (17). Cyr61 is ah extracel- 
lular matrix signaling molecule that promotes cell adhesion, 
proliferation, migration, angiogenesis, and tumor growth (18, 
19). nov (nephroblastoma overexpressed) is an immediate 
early gene associated with quiescence and found altered in 
Wilms tumors (20). The proteins of the CCN family share 
functional, but not sequence, similarity to Wnt-1. All are 
secreted, cysteine-rich heparin binding glycoproteins that as- 
sociate with the cell surface and extracellular matrbc. 

WISP proteins exhibit the modular architecture of the CCN 
family, characterized by four conserved cysteine-rich domains 
(Fig. 3B) (21). The N-terminal domain, which includes the first . 
12 cysteine residues, contains a consensus sequence (GCGC- 
CXXC) conserved in most insulin-like growth factor (IGF)- 
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Fig. 3. (A) Encoded amino acid sequence alignment of human 
WISPs. The cysteine residues of WISP-1 and WISP-2 that are not 
present in WISP-3 are indicated with a dot. (B) Schematic represen- 
tation of the WISP proteins showing the domain structure and cysteine 
residues (vertical lines). The four cysteine residues in the VWC domain 
that are absent in WISP-3 are indicated with a dot. (C) Expression of 
WISP mRNA in human tissues. PCR was performed on human 
multiple-tissue cDNA panels (GLONTECH) from the indicated adult 
and fetal tissues. 

binding proteins (BP). This sequence is conserved in WISP-2 
and WISP-3, whereas WISP-1 has a glutamine in the third 
position instead of a glycine. CTGF recently has been shown 
to specifically bind IGF (22) and a truncated nov protein 
lacking the IGF-BP domain is oncogenic (23). The von Wil- 
lebrand factor type C module (VWC), also found in certain 
collagens and mucins, covers the next 10 cysteine residues, and 
is thought to participate in protein complex formation and 
oligomerization (24). The VWC domain of WISP-3 differs 
from all CCN family members described previously, in that it 
contains only six of the 10 cysteine residues (Fig. 3 A and B). 
A short variable region follows the VWC domain. The third 
module, the thrombospondin (TSP) domain is involved in 
binding to sulfated glycoconjugates and contains six cysteine 
residues and a conserved WSxCSxxCG motif first identified in 
thrombospondin (25). The C-terminal (CT) module contain- 
ing the remaining 10 cysteines is thought to be involved in 
dimerization and receptor binding (26). The CT domain is 
present in all CCN family members described to date but is 
absent in WISP-2 (Fig. 3 A and B). The existence of a putative 
signal sequence and the absence of a transmembrane domain 
suggest that WISPs are secreted proteins, an observation 
supported by an analysis of their expression and secretion from 
mammalian cell and baculovirus cultures (data not shown). 

Expression of WISP mRNA in Human Tissues. Tissue- 
specific expression of human WISPs was characterized by PCR 
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analysis on adult and fetal multiple tissue cDNA panels. 
WlSP'l expression was seen in the adult heart, kidney, lung, 
pancreas, placenta, ovary, small intestine, and spleen (Fig. 3C). 
Little or no expression was detected in the brain, liver, skeletal 
muscle, colon, peripheral blood leukocytes, prostate, testis, or 
thymus. WISP-2 had a more restricted tissue expression and 
was detected in adult skeletal muscle, colon, ovary, and fetal 
lung. Predominant expression of WISPS was seen in adult 
kidney and testis and fetal kidney. Lower levels of WISP-3 
expression were detected in placenta, ovary, prostate, and 
small intestine. 

In Situ Localization of WISP-l and WISP-l. Expression of 
WISP-1 and WISP-2 was assessed by in situ hybridization in 
mammary tumors from Wnt-1 transgenic mice. Strong expres- 
sion of WISP-1 was observed in stromal fibroblasts lying within 
the fibrovascular tumor stroma (Fig. 4 ^-D). However, low- 
level WISP-1 expression also was observed focally within tumor 
cells (data not shown). No expression was observed in normal 
breast. Like WISP-lj WISP-2 expression also was seen in the 
tumor stroma in breast tumors from Wnt-1 transgenic animals 
(Fig. 4 E-H), However, WISP-2 expression in the stroma was 
in spindle-shaped cells adjacent to capillary vessels, whereas 




Fig. 4. {A , C, £, and G) Representative hematoxylin/eosin-stained 
images from breast tumors in Wnt-1 transgenic mice. The correspond- 
ing dark-field images showing WISP-I expression are shown in B and 
/). The tumor is a moderately well-differentiated adenocarcinoma 
showing evidence of adenoid cystic change. At low power (A and B), 
expression of WiSP-I is seen in the delicate branching fibrovascular 
tumor stroma (arrowhead). At higher magnification, expression is seen 
in the stromal(s) fibroblasts (C and £>), and tumor cells are negative. 
Focal expression of WISP-1 ^ however, was observed in tumor cells in 
some areas. Images of WISP-2 expression are shown in E-H, At low 
power (£ and F), expression of WISP-2 is seen in cells lying within the 
fibrovascular tumor stroma. At higher magnification, these cells 
appeared to be adjacent to capillary vessels whereas tumor cells are 
negative (G and H). 



the predominant cell type expressing WISP-1 was the stromal 
fibroblasts. 

Chromosome Localization of the WISP Genes. The chro- 
mosomal location of the human WISP genes was determined 
by radiation hybrid mapping panels. WISP-I is approximately 
3.48 cR from the meiotic marker AFM259xc5 [logarithm of 
odds (lod) score 16.31] on chromosome 8q24.1 to 8q24.3, in the 
same region as the human locus of the novH family member 
(27) and roughly 4 Mbs distal to c-myc (28). Preliminary fine 
mapping indicates that WISP-1 is located near D8S1712 STS. 
WISP-2 is linked to the marker SHGC-33922 (lod = 1,000) on 
chromosome 20ql2-20ql3.1. Human WISP-3 mapped to chro- 
mosome 6q22-6q23 and is linked to the marker AFM211ze5 
(lod = 1,000). WISP-3 is approximately 18 Mbs proximal to 
CTGF and 23 Mbs proximal to the human cellular oncogene 
MYB (27, 29). 

Amplification and Aberrant Expression of WISPs in Human 
Colon Tumors. Amplification of protooncogenes is seen in 
many human tumors and has etiological and prognostic sig- 
nificance. For example, in a variety of tumor types, c-myc 
amplification has been associated with malignant progression 
and poor prognosis (30). Because WISP-I resides in the same 
general chromosomal location (8q24) as c-myc, we asked 
whether it was a target of gene amplification, and, if so, 
whether this amplification was independent of the c-myc lociis. 
Genomic DNA from human colon cancer cell lines was 
assessed by quantitative PGR and Southern blot analysis. (Fig. 
5 A and B). Both methods detected similar degrees of WISP-I 
amplification. Most cell lines showed significant (2t to 4-fold) 
amplification, with the HT-29 and WiDr cell lines demonstrat- 
ing an 8-fold increase. Significantly, the pattern of amplifica- 
tion observed did hot correlate with that observed for c-myc, 
indicating that the c-myc gene is not part of the amplicon that 
involves the WISP-I locus. 

We next examined whether the WISP genes were amplified 
in a panel of 25 primary human colon adenocarcinomas. The 
relative WISP gene copy number in each colon tumor DNA 
was compared with pooled normal DNA from 10 donors by 
quantitative PGR (Fig. 6). The copy number of WISP-1 and 
WISP-2 was significantly greater than one, approximately 
2-fold for WISP-1 in about 60% of the tumors and 2- to 4-fold 
for WISP'2 in 92% of the tumors [P < 0.001 for each). The 
copy number for WISP-3 was indistinguishable from one {P - 
0.166). In addition, the copy number of WISP-2 was signifi- 
cantly higher than that of WISP-1 (P < 0,001). 

The levels of WISP transcripts in RNA isolated from 19 
adenocarcinomas and their matched normal mucosa were 




Fig. 5. Amplification of WISP-1 genomic DNA in colon cancer cell 
lines, (A) Amplification in cell line DNA was determined by quanti- 
tative PCR. (B) Southern blots containing genomic DNA (10 /ig) 
digested with EcoRl (WISP-1) oiXbal (c-myc) were hybridized with 
a 100-bp human WISP-l probe (amino acids 186-219) or a human 
c-myc probe (located at bp 1901-2000). The WISP and myc genes are 
detected in normal human genomic DNA after a longer film exposure. 
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Fig. 6. Genomic amplification of WISP genes in human colon 
tumors. The relative gene copy number of the WISP genes in 25 
adenocarcinomas was assayed by quantitative PCR, by comparing 
DNA from primary human tumors with pooled DNA from 10 healthy 
donors. The data are means ± SEM from one experiment done in 
triplicate. The experiment was repeated at least three times. 

assessed by quantitative PCR (Fig. 7). The level of mSP-l 
RNA present in tumor tissue varied but was significantly 
increased (2- to >25-foId) in 84% (16/19) of the human colon 
tumors examined compared with normal adjacent mucosa. 
Four of 19 tumors showed greater than 10-fold overexpression. 
In contrast, in 79% (15/19) of the tumors examined, WISP-2 
RNA expression was significantly lower in the tumor than the 
mucosa. Similar to WISP-1, WISPS RNA was overexpressed in 
63% (12/19) of the colon tumors compared with the normal 
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Fig. 7, WISP RNA expression in primary human colon tumors 
relative to expression in normal mucosa from the same patient. 
Expression of WISP mRNA in 19 adenocarcinomas was assayed by 
quantitative PCR. The Dukes stage of the tumor is listed under the 
sample number. The data are means ± SEM from one experiment 
done in triplicate. The experiment was repeated at least twice. 



mucosa. The amount of overexpression of WISP-3 ranged from 
4- to >40-fold. 

DISCUSSION 

One approach to understanding the molecular basis of cancer 
is to identify differences in gene expression between cancer 
cells and normal cells. Strategies based on assumptions that 
steady-state mRNA levels will differ between normal and 
malignant cells have been used to clone differentially ex- 
pressed genes (31). We have used a PCR-based selection 
strategy, SSH, to identify genes selectively expressed in 
C57MG mouse mammary epithelial cells transformed by 
Wnt-1. 

Three of the genes isolated, WISP-l, WISP-2, and WISP-3, 
are members of the CCN family of growth factors, which 
includes CTGF, Cyr61, and nov, a family not previously linked 
to Wnt signaling. 

Two independent experimental systems demonstrated that 
WISP induction was associated with the expression of Wnt-1. 
The first was C57MG cells infected with a Wnt-1 retroviral 
vector or C57MG cells expressing Wnt-1 under the control of 
a tetracyline-repressible promoter, and the second was in 
Wnt-1 transgenic mice, where breast tissue expresses Wnt-1, 
whereas normal breast tissue does not. No WISP RNA expres- 
sion was detected in mammary tumors induced by polyoma 
virus middle T antigen (data not shown). These data suggest 
a link between Wnt-1 and WISPs in that in these two situations, 
WISP induction was correlated with Wnt-1 expression. 

It is not clear whether the WISPs are directly or indirectly 
induced by the downstream components of the Wnt-1 signaling 
pathway (i.e., ^-catenin-TCF-l/Lefl). The increased levels of 
mSP RNA were measured in Wnt-l-transformed cells, hours 
or days after Wnt-1 transformation. Thus, WISP expression 
could result from Wnt-1 signaling directly through p-catenin 
transcription factor regulation or alternatively through Wnt-1 
signaling turning on a transcription factor, which in turn 
regulates WISPs. 

The WISPs define an additional subfamily of the CCN family 
of growth factors. One striking difference observed in the 
protein sequence of WISP-2 is the absence of a CT domain, 
which is present in CTGF, Cyr61, nov, WISP-1. and WISP-3. 
This domain is thought to be involved in receptor binding and 
dimerization. Growth factors, such as TGF-/3, platelet-derived 
growth factor, and nerve growth factor, which contain a cystine 
,knot motif exist as dimers (32). It is tempting to speculate that 
WISP-1 and WISP-3 may exist as dimers, whereas WISP-2 
exists as a monomer. If the CT domain is also important for 
receptor binding, WISP-2 may bind its receptor through a 
different region of the molecule than the other CCN family 
members. No specific receptors have been identified for CTGF 
or nov. A recent report has shown that integrin avjSa serves as 
an adhesion receptor for Cyr61 (33). 

The strong expression of WISP- 1 and WISP-2 in cells lying 
within the fibrovascular tumor stroma in breast tumors from 
Wnt-1 transgenic animals is consistent with previous obser- 
vations that transcripts for the related CTGF gene are pri- 
marily expressed in the fibrous stroma of mammary tumors 
(34). Epithelial cells are thought to control the proliferation of 
connective tissue stroma in mammary tumors by a cascade of 
growth factor signals similar to that controlling connective 
tissue formation during wound repair. It has been proposed 
that mammary tumor cells or inflammatory cells at the tumor 
interstitial interface secrete TGF-j31, which is the stimulus for 
stromal proliferation (34). TGF-)31 is secreted by a large 
percentage of malignant breast tumors and may be one of the 
growth factors that stimulates the production of CTGF and 
WISPs in the stroma. 

It was of interest that WISP-1 and WISP-2 expression was 
observed in the stromal cells that surrounded the tumor cells 
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(epithelial cells) in the Wnt-l transgenic mouse sections of 
breast tissue. This finding suggests that paracrine signaling 
could occur in which the stromal cells could supply WISP-1 and 
WISP-2 to regulate tumor cell growth on the WISP extracel- 
lular matrix. Stromal cell-derived factors in the extracellular 
matrix have been postulated to play a role in tumor cell 
migration and proliferation (35). The localization of WISP-l 
and WISP-2 in the stromal cells of breast tumors supports this 
paracrine model. 

An analysis of WISP-] gene amplification and expression in 
human colon tumors showed a correlation between DNA 
amplification and overexpression, whereas overexpression of 
WlSP-3 RNA was seen in the absence of DNA amplification. 
In contrast, WISP-2 DNA was amplified in the colon tumors, 
but its mRNA expression was significantly reduced in the 
majority of tumors compared with the expression in normal 
colonic mucosa from the same patient. The gene for human 
WISP-2 was localized to chromosome 20ql2-20ql3, at a region 
frequently amplified and associated with poor prognosis in 
node negative breast cancer and many colon cancers, suggest- 
ing the existence of one or more oncogenes at this locus 
(36-38). Because the center of the 20ql3 amplicon has not yet 
been identified, it is possible that the apparent amplification 
observed for WISP-2 may be caused by another gene in this 
amplicon. 

A recent manuscript on rCop-I, the rat orthologue of 
WISP-2, describes the loss of expression of this gene after cell 
transformation, suggesting it may be a negative regulator of 
growth in cell lines (16). Although the mechanism by which 
WISP-2 RNA expression is down-regulated during malignant 
transformation is unknown, the reduced expression of WISP-2 
in colon tumors and cell lines suggests that it may function as 
a tumor suppressor. These results show that the WISP genes 
are aberrantly expressed in colon cancer and suggest that their 
altered expression may confer selective growth advantage to 
the tumor. 

Members of the Wnt signaling pathway have been impli- 
cated in the pathogenesis of colon cancer, breast cancer, and 
melanoma, including the tumor suppressor gene adenomatous 
polyposis coli and ^-cateriin (39). Mutations in specific regions 
of either gene can cause the stabilization and accumulation of 
cytoplasmic )3-catenin, which presumably contributes to hu- 
man carcinogenesis through the activation of target genes such 
as the WISPs, Although the mechanism by which Wnt-l 
transforms cells and induces tumorigenesis is unknown, the 
identification of WISP^ as genes that may be regulated down- 
stream of Wnt-l in C57MG cells suggests they could be 
important mediators of Wnt-l transformation. The amplifica- 
tion and altered expression patterns of the WISP^ in human 
colon tumors may indicate an important role for these genes 
in tumor development. 
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